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ABSTRACT 

We investigate the biases and uncertainties in estimates of physical parame- 
ters of high-redshift Lyman break galaxies (LBGs), such as stellar mass, mean 
stellar population age, and star formation rate (SFR), obtained from broad-band 
photometry. These biases arise from the simplifying assumptions often used in 
fitting the spectral energy distributions (SEDs). By combining ACDM hierar- 
chical structure formation theory, semi-analytic treatments of baryonic physics, 
and stellar population synthesis models, we construct model galaxy catalogs from 
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which we select LBGs at redshifts z ~ 3.4, 4.0, and 5.0. The broad-band photo- 
metric SEDs of these model LBGs are then analysed by fitting galaxy template 
SEDs derived from stellar population synthesis models with smoothly declining 
SFRs. We compare the statistical properties of LBGs' physical parameters - such 
as stellar mass, SFR, and stellar population age - as derived from the best-fit 
galaxy templates with the intrinsic values from the semi-analytic model. We find 
some trends in these distributions: first, when the redshift is known, SED- fitting 
methods reproduce the input distributions of LBGs' stellar masses relatively well, 
with a minor tendency to underestimate the masses overall, but with substantial 
scatter. Second, there are large systematic biases in the distributions of best-fit 
SFRs and mean ages, in the sense that single-component SED-fitting methods 
underestimate SFRs and overestimate ages. We attribute these trends to the 
different star formation histories predicted by the semi-analytic models and as- 
sumed in the galaxy templates used in SED-fitting procedure, and to the fact that 
light from the current generation of star-formation can hide older generations of 
stars. These biases, which arise from the SED-fitting procedure, can significantly 
affect inferences about galaxy evolution from broadband photometry. 

Subject headings: galaxies: evolution — galaxies: fundamental parameters — 
galaxies: high-redshift — galaxies: statistics — galaxies: stellar content — meth- 
ods: statistical 



1. Introduction 



Over the last decade, astronomy has experienced a boom of panchromatic surveys of 
high redshift galaxies - such as the Great Observatories Origins Dee p Survey (GOODS , 
Giavalisco et al.ll2004al ) and the Cosmic Evolution Survey (COSMOS. IScoville et al.l 120071 ) 
- thanks to the developments of space telescopes, like Hubble Space Telescope (HST) and 
Spitzer Space Telescope (Spitzer), and also of ground-based facilities, including Keck obser- 
vatory. Very Large Telescope (VLT), and Subaru telescope. Spectral energy distributions 
(SEDs) constructed from the photometric data of wide wavelength coverage have been used 
to constrain galaxy parameters, such as stellar masses, SFRs, and stellar population ages, 
of hi gh redshift galaxies yia comparison with simple stellar population synthesis models 



Papovich et all l200ll : IShaplev et all l200ll . l2005l : IShim et al.l l2007l : IStark et all 12007 



Verma et al]l2007h . Accompanied byselection techniques to separate high-redshift galaxy 
samples efficiently via their broad-band colors and magnitudes, these surveys have provided 
increasing information about the nature of high-redshift galaxies. 
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However, despite recent advances in our knowledge of the properties of high redshift 
galaxies, more accurate estimation of physical parameters is necessary for addressing sev- 



eral important issues in galaxy evolution and cosn iology. For exam 



certain whether the Lyman break galaxies (LBGs, ISteidel et al.l Il996l : iGiavaliscd |2002| ) - 



j le, it is as yet u n- 



high-redshift, star-forming galaxies selec ted according to their rest-frame ultravio let (UV) 
colors - will evolve into large ellipticals (lAdelberger et al.lll998l : ISteidel et al.lll998l ). or into 
smaller spheroids, such as small ellipticals /spiral bulges or subgalactic s tructures, which will 
later be merged to form larger galaxies at 2; ~ (ILowenthal et al.lll997l : ISawicki fc Yedll998 



Somerville et al.ll200ll ). Accurate knowledge of physical parameters, such as stellar masses 
and SFRs, serves as one of the important elements in discriminating among the possible evo- 
lutionary descendants of LBGs. Also, accurately constraining LBG's stellar-population ages 
and stellar masses is necessary for determining how much the LBG populations contributed 
to cosmic reionization. Better estimation of these parameters is also crucial in comparing 
different galaxy populations at different redshifts. 

The stellar populations of LBGs as well as other populations of high-redshift galaxies 
have been studied by several authors using SED-fitting methods to compare the photometric 
SEDs of observed galaxi e s wit h various galaxy spectra from stellar population synthesis 



models. iPapovich et al.l (I2OOII ) investigated 33 spectroscopically-confirmed LBGs with a 
redshift range, 2.0 < z < 3.5. They found that the mean lower limit of stellar masses 
of these LBGs is ~ 6 X 10^ Mo with upper limits of ~ 3-8 ti r nes la rger and that the 



mean age is ~ 120 Myr (assuming solar metallicity and ISalpeterl (119551 ) IMF) with broad 
range between 30 Myr and 1 Gyr. More importantly, they concluded that the most robust 
parameter constrained through SED-fitting is stellar mass, while stellar population age and 
especially star formation rate (SFR) are only poorly constrained. Also, they speculated that 
SED-fitting methods (single- component fitting) can only give lower limits to galaxies' stellar 
masses because the results from the SED-fitting r nethods are largely dr iven by the light from 
the most massive, most recently formed stars. IShapley et al.l (1200 ll ) studied the physical 
parameters of 74, 2; ~ 3 LBGs with spectroscopic redshifts. They confirmed that among 
various physical parameters, stellar mass is the most tightly constrained, and that constraints 
on other parameters such as dust extinction and age are weak. With the assumption of solar 
metallicity and a Salpeter IMF, the median age (defined as time since the onset of current 
star formation) was found to be t^/ ~ 320 Myr with a large spread from several Myr to 
mo re than 1 Gyr. The m edian stellar mass w as 1.2x10-^° Mq , which is higher than 
the iPapovich et al.l (1200 ll ) value. However, the IShapley et al.l (1200 ll ) LBGs were generally 
brighter than LBGs in the IPapovich et al.l (120011 ) sample. The two studies found si milar 
stellar masses for LBGs with similar rest-frame UV luminosities. IShapley et al.l (120011 ) also 
found that about 20 % of the LBGs have best-fit ages older than 1 Gyr, stellar masses larger 
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than 10^° Mq, and SFR ~ 30 Mq yr~^, and they interpreted these a s galaxies which have 
formed their stars over a relatively long period in a quiescent manner. IShapley et al.l (120011 ) 
also noted that more luminous galaxies are dustier than less luminous ones and that younger 
galaxies are dustier and have higher SFR. 



More recently, IVerma et al.l (120071 ) studied 21 2; ~ 5 LBGs, six of which have confirmed 
spectroscopic redshifts and the remaining 15 of which have photometric redshifts. These 
LBGs were found to be moderately massive with median stellar mass ~ 2 x 10^ Mq, and 
to have high SFRs with a median of ~ 40 MQyr~^. They also found that the stellar mass 
estimates are the most robust of all derived properties. Best-fit ages have a large spread with 
a me dian value ^ 25 My r, assuming a metallicity of one- fifth solar (0.2 Zc?,) and a Salpeter 
IMF. IVerma et al.l (120071 ) also compared their z ~ 5 LBGs with lShapley et al.l (l200ll )'s z ~ 3 



LBGs, assuming the same IMF and metallicity (i.e. solar metallicity and Salpeter IMF), and 
concluded that these two samples of LBGs with similar rest-frame UV luminosity are clearly 
different. Specifically, 2; ~ 5 LBGs are much younger (< 100 Myr) and have lower stellar 
masses (~ 10^ Mq) than z ~ 3 LBGs. The fraction of young (age < 100 Myr) galaxies is 
~ 70% for 2; ~ 5 LBGs while it is ~ 30% for 2; ~ 3 LBGs. They also concluded that these 
2 ~ 5 LBGs are the likely progenitors of the spheroidal co mponents of present -day massive 



galaxies, based on their high stellar mass surface densities. iStark et al.l (120071 ) analysed 72 
z ~ 5 galaxies (not just LBGs) with photometric redshifts. They performed an SED-fitting 
analysis on a spectroscopically confirmed subset of 14 2; ~ 5 galaxies to derive best-fit stellar 
masses ranging between 3 x 10*^ and 2 x 10^^ Mq, and ages from 1 Myr to 1.1 Gyr. Three out 
of these 14 galaxies have stellar masses in excess of 10^^ Mq. Using stellar mass estimates 
from both spectroscopically confirmed and candidate galaxies, they calculated a stellar mass 
density at z ~ 5 of 6 x 10^ Mq Mpc~^, which is much larger than the integration of the 
star formation rate density (SFRD) from z ~ 10 to 5. They attributed this discrepancy 
either to significan t dust extinction or to undetected low-luminosity star-forming galaxies at 
z > 5. IShim et al.l (120071 ) studied the properties of 1088 massive LBGs whose best-fit stellar 
masses are larger than 10^^ Mq at z ~ 3. They derived stellar masses of these massive LBGs 
through SED-fitting with the assumption of a Salpeter IMF, and noted that LBGs which 
are detected in mid-infrared wavelength are the ones with large stellar mass and severe dust 
extinction among LBG population. 

The SED-fitting meth od is also used to constr ain physical parameters for high -z Lyman 



a ern it ting galaxies (LAEs. iFinkelstein et al.ll2007l ) and for submillimetre galaxies (IDye et al. 
20081 ). IFinkelstein et al.l (120071 ) analysed z ~ 4.5 LAEs through SED fitting and found their 



ages range from 1 Myr to 200 Myr (a ssuming a constan t star formation history), and stellar 



masses from 2 x 10^ — 2 x 10^ Mq. iDye et al.l (120081 ) used SED fitting to study physical 



parameters of SCUBA (Submillimetre Common-User Bolometer Array) sources, finding an 
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average stellar mass of 



IQll.SiO.lj^^ 



0- 



However, the estimated parameter values and implications derived from these values 
are prone to several sources of error including errors inherent in SED-fitting methods. It 
is unclear whether or not SED-fitting methods deliver biased estimates of parameters such 
as stellar mass, age, and SFR. It is also unclear how much worse the parameter estimates 
become if spectroscopic redshifts are unavailable and the photometric data must be used to 
constrain not only the star formation history, but also the redshift. The critical question is 
how far we can trust the physical parameters derived using mainly photometric data under 
the assumption of simple star- formation histories (SFHs), when real galaxies are expected 
to have more complex formation/evolution histories with several possible episodes of star- 
formation. 

In this paper, we address this issue by comparing the statistical distributions of intrin- 
sic physical parameters of model Lyman break galaxies (LBGs) from semi-analytic models 
(SAMs) of galaxy formation, with the estimates of those same parameters derived from 
stellar population synthesis models through the commonly used SED-fitting method. 

There have been several studies which tried to constrain the stellar populations of LBGs 



as semi-analytic models (e.g._ 
drodynamic simu 



Somerville et a 


2001 




Idzi et al. 


2004 


or cosmological hy 


Naeamine et al. 


2OO5I; 


Nieht et al. 


20061; 


Finlator et al. 


2007) 



Somerville et al.l (1200 ll ) compared z ~ 3 and z ~ 4 observed LBGs w ith semi-analytic mod- 



els to investigate possible scenarios for LBGs, and lldzi et al.l (120041 ) put constraints on the 
proper ties of z ~ 4 LBGs thro ugh comparison of colo rs of observed LBGs and SAMs model 



LBGs. iNagamine et al.l (120051 ) and lNight et al.l (120061 ) used cosmological simul ations to place 
constraints on properties of UV-selected, 2 ~ 2 galaxies and ^ ~ 4 — 6 LBGs. iFinlator et al. 



( 120071 ) constructed a set of galaxy templates from hydrodynamical simulations and used 
them to place constraints on properties of 6 2 < 5.5 observed galaxies. Our approach is 
distinct from these studies. Here, instead of comparing model galaxies with observed ones, 
we perform SED fitting on model galaxies, derive best-fit parameters, and compare them 
with the intrinsic values given in the model, trying to understand quantitatively as well as 
qualitatively the biases and uncertainties of SED-fitting methods in constraining the physical 
parameters of LBGs. 

To accomplish this, we construct galaxy formation histories by combining ACDM hier- 
archical structure formation theory with semi-analytic treatments of gas cooling, star for- 
mation, supernova feedback, and g alaxy mergers. Then, usin g the stellar population synthe- 
sis models of Bruzual & Chariot (IBruzual &: CharlotI l2003l . hereafter BC03), and a simple 
model for dust extinction, we derive the photometric properties of these model galaxies. 
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After selecting LBGs through appropriate color criteria, we compare the photometric spec- 
tral energy distributions (SEDs) in observed-frame optical (HST/ACS) through mid-infrared 
[Spitzer /YRAG) passbands of these model LBGs with those of SED templates from stellar 
population synthesis model of BC03 using a x^-Kiinimization method. In this SED-fitting 
procedure, a large multi-parameter space is explored to minimize any prior. Each parameter 
- such as star-formation time scale, age, dust-extinction, and metallicity - in this parame- 
ter space spans a broad as well as physically realistic range. The distributions of physical 
parameters of these model LBGs derived from this SED-fitting method are compared with 
the input distributions of these parameters from the semi-analytic model. 

Various parameters in the semi-analytic model are calibrated to reproduce the col- 
ors of observed LBGs in the southern field of the GOODS (GOODS-S). The GOODS 



( iGiavalisco et al.l l2004al ) is a deep, multiwavelength survey which covers a combined area 
of ~ 320 arcmin^ in two fields - GOODS-N centered on the Hubble Deep Field-North 
(HDF-N), and GOODS-S centered on the Chandra Deep Field-South (CDF-S). The exten- 
sive wavelength coverage of the photometric data in the GOODS fields including HST/ACS 
[Advanced Camera for Surveys) and Spitzer /IUAC {Infrared Array Camera) is beneficial 
for the derivation of various physical parameters of galaxies including stellar mass, SFR, 
and age. Deep observations (reaching z^^q ~ 26.7) in the GOODS fields can probe the 
high redshift universe in a comprehensive and statistically meaningful manner. This has en- 



Giava. 



2006 



isco et all 2004b 



Yan et al.l 



20061 : 



Idzi et ahlbooi ] 


^apovich et al. 


2004 


Lee et al.ll2006: 


Ravindranath et al. 


Stark et al. 


200' 


1\ 


Verma et al. 


2007 


). In this work, we use band- 



passes for the HST/ACS fihers F435W (S435), F606W (V^eoe), F775W (2775), F850LP (zgso), 
VLT/ISAAC (Infrared Spectrometer And Array Camera) J, H, Ks bands, and Spitzer /YRAC 
3.6, 4.5, 5.8, 8.0 /im channels for SED-fitting. The CTIO (Cerro Tololo Inter-American Ob- 
servatory) MOSAIC U band is used for selecting U-dropouts. 

A description of the semi-analytic model of galaxy formation used in this work is given in 
§ 2, with the Lyman break galaxy sample selection and the SED-fitting procedures explained 
in § 3. The statistical properties of the derived physical parameters for model LBGs are 
shown in § 4. We analyse, in detail, the effects of various factors on the biases in SED-fitting 
in § 5, and the bias in estimating SFR from rest-frame UV magnitude in § 6. We discuss 
the effects of these biases on the galaxy evolution studies in § 7 and a summary and our 
conclusions are given in § 8. Throughout the paper, we adopt a fiat ACDM cosmology, with 
(f2m,^A) = (0.3,0.7), and Hq = lOO/i kr n Mpc ~^, where h=0.7, and all magnitudes are 
expressed in the AB magnitude system (IOkelll974j ). 
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2. Semi-Analytic Models of Galaxy Formation 



Semi-analytic models (SAMs) of galaxy formation are embedded within the framework 
of a ACDM initial power spectrum and the theory of the growth and collapse of fluctuations 
through gravitational instability. The models include simplifled physical treatments of gas 
cooling, star formation, supernova feedback, dust extinction, and galaxy merging. In this 
work, we use the SAMs model ru n which was constructed and used in lldzil ( 120071 ). The 
Idzj (|2007[) work was base d on the ISomerville &: Primackl (119991 ) model, which is also used 
bv ISomerville et al.l (120011 ) and lldzi et al.l (120041 ). Speciflcally, among various model runs in 
Idzil ( 120071 ). we use here the model run which showed the closest match to the rest-frame 
UV-continuum and UV-optical color^ of observed U- & B-dropouts in the GOODS-S fleld. 
The details of the model can be found in the above references, and here we briefly highlight 
the most important aspects of the model (including some changes from the models described 
in the original papers). 



2.1. Description of Semi- Analytic Models 



The semi-analytic model used in this work implements the method of lSomerville fc Kolatt 



(119991 ) to build merging histories of dark matter halos. After dark matter merging histo- 
ries are constructed, the semi-analytic treatments for various physical processes, such as gas 
cooling, mergers, star formation (both in a merger-induced burst mode and in a quiescent 
mode), supernovae feedback, chemical evolution, and dust extinction, are applied to realize 
predictions of the formation history of a statistical ensemble of galaxies. 

A newly formed dark matter halo (residing at top of the tree) contains pristine shock- 
heated hot gas at the virial temperature. When a halo collapses or undergoes a merger with 
a larger halo, the associated gas is assumed to be shock-heated to the virial temperature 
of the halo. This gas then radiates energy and consequently cools. For small halos at high 
redshift, cooling is limited by the accretion rate, since the amount of gas that can cool at any 
given time cannot exceed the amount of hot gas contained within the halo's virial radius. 

Once halos merge, the galaxies within them remain distinct for some time. The central 
galaxy of the largest progenitor halo is set as the central galaxy of the merged dark matter 
halo. All the other galaxies become sa tellites, which then fall in towards the central galaxy 



due to dynamical friction. Unlike in ISomerville et al.l (1200 ll ). collisions between satellite 



^^606 - *775 and Vfjoe - IRAC 3.6 /^m colors are used for U-dropouts, and 1775 - zgso and 2^50 - IRAC 4.5 
/im colors are used for B-dropouts as diagnostics to determine best-fit model run. 
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galaxies are neglected here. 

Both quiescent and merger-driven modes of star formation are included in the model. 
We assume that every "major" galaxy-galaxy merger above a certain mass-ratio threshold 
(1:4) triggers a starburst, which converts 100 % of the available gas into stars in a tenth of 
a halo dynamical time. Quiescent star formation is modeled with a Kennicutt-like law, such 
that the star formation rate is proportional to the mass of cold gas in the disk divided by 
the dynamical time. We also scale the SFR by a power-law function of the galaxy circular 
velocity, such that star formation is less efficient in low-mass galaxies. This mimics the effect 



of a S F threshold, as implemented in more recent semi-analytic models (e.g. ISomerville et al. 



20081 ). The total star formation rate of a galaxy is the sum of the burst and quiescent modes. 



Supernova feedback is also modeled via a simple recipe, where the rate of reheating of 
cool gas by supernovae is proportional to the SFR times a power-law function of circular 
velocity, such that the reheating is more efficient in smaller mass galaxies. If the halo's virial 
velocity is less than a preset ejection threshold (set to 100 km in our model), then all of 
the reheated gas is ejected from the halo; otherwise the reheated gas is placed in the hot gas 
reservoir within the halo. The gas and metals that are ejected from the halo are distributed 
outside of the halo with a continuation of the isothermal profile that we assumed inside 
the halo. This material falls in gradually as the virial radius of the halo increases due to the 
falling background density of the Universe. 

Chemical evolution is treated by assuming a constant mean mass of metals produced 
per mass of stars. The metals produced in and ejected from stars are first deposited into the 
surrounding cold gas, at which point they may be ejected from the disk and mixed with the 
hot halo gas. The metallicity of any newly formed stars is set to equal the metallicity of the 
ambient cold gas at the time of formation. 

For each galaxy, the SAM predicts the two-dimensional distribution of the mass in 
stars of a given age and metallicity. We convolve this distribution with the SSP models of 



BC03 (using the Padova 1994 isochrones) to create a synthetic SED. We adopted a lChabrier 



( 120031 ) IMF with lower ar id uppe r mass cutoffs of = 0.1 Mq and mu = 100 M( 



0- 

We then use the model of iMadaul (119951 ) to account for the opacity of intervening HI in the 
intergalactic medium as a function of redshift, and convolve the SEDs with the response 
functions appropriate to the ACS, ISAAC, IRAC, and MOSAIC photometric bands. 

Dust extinction is treated assuming that the face-on optical depth in the V-band is 
Tvfl = Tdustfi X (m*)^'^"°* where Tdust,o and Pdust are free parameters set as 1.2 and 0.3, to match 
the observatio n s in t he GOODS-S. This choice is motivated by the observational results of 



Hopkins et al.l (120011 ). The dependence of the extinction on wavelength (the attenuation 
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curve) is calculated using a Calzetti attenuation curve (ICalzetti et al.ll2000l ). 

In summary, we use the semi-analytic models to build galaxies that have rich and varied 
star formation histories that are motivated by the hierarchical galaxy formation picture. Star 
formation is bursty and episodic. Stars have a distribution of metals consistent with the 
assumed star formation history, and the dust content is plausible. In contrast, the simple 
SED templates which we will use to analyze these model galaxies have uniform metallicity 
and monotonically declining star formation rates. 



3. Lyman Break Galaxy Samples and Spectral Energy Distribution Fitting 
3.1. Model Galaxy Catalog and Lyman Break Color Selection 



We created a suite of model runs by varying the uncertain model parameters controlling 
the burstiness of star formation and dust extinction. As stated in § [2], the rest-frame UV- 
continuum and UV-optical colors of model LBGs from each semi-analytic model run were 
compared with the colors of LBGs observed in the GOODS-S field, and the best-fit model 
run was selected based on these comparisons. From this best-fit model run, a model galaxy 
catalog, which contains a total of 44281 galaxies within a redshift range 2.3 < z < 5.7, has 
been constructed. This catalog carries the various physical parameters such as stellar mass 
and SFR, along with the (dust-extinguished) broad band photometry for several bandpasses 
used in the GOODS observations. The photometric data in this catalog are derived by 
combining the SAM's star formation and enrichment histories with BC03 stellar population 
synthesis models, and including models for dust extinction and absorption by the IGM, as 
described in § [2l 

These photometric data are used to select high-redshift, star-forming galaxies through 
the Lyman-break color selection technique. The Lyman-break color selection technique - 
which uses the 'Lyman-break' feature at A ~ 912 A in galaxy spectra and the Lyman-a 
forest flux deficit between 912 A and 1216 A together with blue colors at longer rest-frame 
UV wavelengths to identify star-forming galaxies located at high redshift - has been shown 
to be an effective way to const ruct large sarnples of high-redshift, sta r-forming galaxies from 
optical photometric data sets (IMadau et al.lll996l : ISteidel et al.ll2003l ). Spectroscopic follow- 
ups for the LBG samples have ver ified the robustness and efficiency of this te chnique in 
building up high- 2; galaxy samples (ISteidel et al.lll996l . Il999l : IVanzella et al.ll2006l ) 



The color criteria used in this work to select U-, B-, & V-dropout model LBGs are as 
follows: 
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{U - 5435) > 0.62 + 0.68 X (5435 - zs5o) A (1) 

(U - 5435) > 1.25 A (2) 

(S435 - ^85o) < 1-93 (3) 

for U-dropouts (z~3.4), 

(E435 - Veoe) > 1-1 + (^606 - 285o) A (4) 

(^435 - ^606) > 1.1 A (5) 

(Veoe - Z85o) < 1.6 (6) 

for B-dropouts (z~4), and 



((^606 - ^775) > 1.4667 + 0.8889 x (2775 - 2:850)) V ((l^eoe - ^775) > 2.0) A (7) 

(Veoe - ^775) > 1.2 A (8) 

(^775 - Z85o) < 1.3 (9) 

for V-dropouts (z~5). 

Here, A means logical 'AND', and V is logical 'OR'. 

The criteria for B- & V-dropouts used here are the same ones which have been tested 
and used in previous works on GOOD S survey fields (e.g.. iGiavalisco et al.ll2004bl : lldzi et al. 
2OO4J : iPapovich et al.ll2004l : lldzil 120071 ). The criteria for U-dropouts are set (1) to select z ~ 
3 star-forming galaxies, (2) to effectively avoid contamination from foreground stars, and 
(3) to minimize the overlap with B-dropouts. These color criteria and additional magnitude 
cuts for the ACS z-band {zab < 26.6) and IRAC channel 1 or channel 2 (m3.6/im < 26.1 or 
"^4.5/im < 25.6), which correspond to detection magnitude limits for each bands in GOODS- 
S, give samples of LBGs containing 2729 objects for U-dropouts with redshift range 3.1 < 
z < 3.6, 2638 objects for B435-dropouts with 3.6 < z < 4.6, and 808 for Veoe-dropouts with 
4.6 < ;z < 5.6. 

Figure [T] shows the redshift distributions of model U-, B-, and V-dropout galaxies se- 
lected with the above criteria. The color selection criteria effectively isolate three distinct 
redshift intervals. Because the models include a realistic distribution of galaxy luminosities, 
the galaxies preferentially lie toward the lower redshift boundary of each redshift interval, 
as they do for real surveys. 
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3.2. Special Energy Distribution Fitting 

The goal of this work is to compare the statistical distributions of various physical 
parameters of model Lyman break galaxies with the ones derived from the SED-fitting 
method. To do this, we performed x^-minimization to find the best-fit stellar population 
synthesis model templates and best-fit parameters derive d from them. H ere, we use as the 
stellar population synthesis model the BC03 model with IChabrien (120031 ) IMF and Padova 
1994 evolutionary tracks. This is the same set of SSP models used to construct the SEDs 
for the SAM galaxies. Only few discrete values of metallicities are available in the BC03 
model library; among these, two sub-solar (0.2 Zq and 0.4 Zq) and solar metallicities were 
used. Studies of high redshift LBG spectra have shown that their rnetallicities are sub solar 
or solar (Z ~ 0.2 - 1.0 Zq) JPettini et ahlboOll : Ishaplev et allbooj: lAndo et al.ll2004h . For 



internal dust extinction, the Calzetti extinction law (iCalzetti et al 



2000) is used with 0.000 



< E{B — V) < 0.950 with a step size 0.025, and intergalactic absorption by neutral hydrogen 
according to iMadaul (119951 ) is included. The star formation histories are parameterized as 
cc e^*/"^ with the e-folding time, r ranging between 0.2 and 15.0 Gyr, and the time since onset 
of star formation, t spanning from 10 Myr to 2.3 Gyr, with limiting t being smaller than 
the age of the universe at each corresponding redshift. The parameter values used in our 
SED-fitting are summarized in table. [H From the observed LBGs in GOODS-S, mean errors 
for different magnitude bins at each passband of ACS, ISAAC and IRAC are calculated, 
and assigned to each SAM galaxy photometric value according to their magnitudes in the 
calculation of x^- Fluxes of the SAM galaxies in the ACS -B435, Vqqq, in^, zs5q bands (-B435 
band photometry are not used in fitting for Veoe-dropouts), ISAAC J, H, Ks bands, and 
IRAC 3.6, 4.5, 5.8, 8.0 /im channels are compared with those of BC03 model template^, 
and based on the calculated values, the best-fit galaxy template is determined for each 
model LBG. The best-fit stellar mass is calculated by multipling the mass-to-light ratio and 
bolometric luminosity of the best-fit template. 

The best-fit stellar mass equals to the integration of the SFR over time, correcting for 
stellar mass recycling. 



Ax [[\{t',T)dt' -r{to,T)] 
Jo 



(10) 



Here, to and r are the best-fit t and r respectively, is the best-fit stellar mass, and 
A is a normalization factor. \E'(t,r) and r[t,T) are unnormalized SFR and stellar mass 



^U-band data is used only for selecting U-dropout samples, and is not used in the SED-fitting procedure. 
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recycling fraction respectively, and can be known from the best-fit t and r. The best-fit 
instantaneous SFR is then, 



where ^'(to,^) = e"*"/^. 

Not correcting for recycling would result in a ~ 29 % underestimation of the correct SFR 
and 8 ~ 9 % overestimation of the mean ag^. The statistical behaviors of various physical 
parameters, such as stellar mass, SFR, and mean age, derived from these best-fit templates 
are compared with the intrinsic distributions from the SAMs in the following sections. 

4. Constraining Physical Parameters using the SED-Fitting Method 



In the ideal case, in which the SED-fitting methods can determine various physical 
parameters accurately, the one-to-one comparison between intrinsic SAM parameters and 
best-fit parameters from SED-fitting should show no deviation from a straight line with a 
slope of unity. However, there are several factors that introduce errors in the parameters 
derived from SED-fitting, which include the complex star formation/merging histories of 
individual galaxies and photometric errors. So, the interesting question is how well the SED- 
fitting methods can recover the intrinsic or true statistical distributions of various physical 
parameters - such as stellar mass, star formation rate (SFR), mean stellar population age, 
and so on - which are used to constrain galaxies' formation histories. 

The star formation rate of a galaxy is an ill-defined quantity. Strictly speaking, it is 
always zero except at those instants in time when a collapsing ball of gas begins to generate 
energy through nuclear fusion. It is common in the literature to implicitly average over some 
time interval At without specifying At. In our view, it is important to be specific about the 
timespan for this averaging. Henceforth in this paper, we define the SFR to be the mass in 
stars formed in the past 100 Myr. Our reason for this choice is as follows: unlike the BC03 
model galaxy templates whose SFRs are assumed to decrease exponentially with time, the 
star formation activity in semi-analytic model galaxies - and probably in real LBGs also 



•^In some studies, SFR is derived first from the best-fit t and r, then stellar mass is calculated by integrating 
the SFR over time. In these cases, failure to correct for stellar mass recycling results in an overestimation 
of stellar mass. 





4.1. Physical Parameters of Model Lyman Break Galaxies 
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e.g. iPapovich et al.ll2005l ) - is more episodic and complex. So, instantaneous SFRs or SFRs 



averaged over very small At are less meaningful for SAM LBGs. 

In this work, we repeat our SED-fitting experiment with several different choices of 
conditions: (1) holding redshift fixed and fitting redshift as an additional free parameter, (2) 
using various combinations of passbands of ACS, ISAAC, and IRAC, and (3) with limited 
T ranges. We also try two-component templates as well as single-component templates. 

In this section, we show the SED-fitting results in cases where we use all the passbands 
from ACS, ISAAC, and IRAC: first, assuming the redshift is known, and second, fitting 
redshift as an additional free parameter. 

Figures U\ El and S] compare the intrinsic distributions to the fitted values of stellar 
mass, SFR (averaged over the past 100 Myr), and (mass- weighted) mean stellar population 
age for z ~ 3.4 (U-dropouts), 4.0 (B-dropouts), and 5.0 (V-dropouts) model LBGs, with the 
redshift fixed to the actual value in the SAM model catalog. 

Figure [2] shows that the SED-fitting method recovers relatively well the input distribu- 
tions of stellar masses in spite of the very different star formation histories (SFHs) in the 
semi-analytic model and in the templates from the BC03 stellar population synthesis models. 
The mean values of the SED-derived stellar masses of U-, B-, and V-dropouts differ from the 
intrinsic mean values by ~ 19 % (U-dropouts), 25 % (B-dropouts), and 25 % (V-dropouts) of 
the intrinsic mean values, in the sense that the SED-fitting method systematically underesti- 
mates the mean values of stellar masses for all of three dropout samples. This trend may be 
attributed to the fact that light from a recent generation of star formation can easily mask 
(some portion of) the presence of an older stellar population. Our results are consistent with 
the earlier arguments that the stellar masses o f LBGs are underestimated when derived f rom 



(single-component) SED-fitting methods (e.g. iPapovich et al.ll200ll : IShapley et al.ll2005l ). 



However, figures [3] and H] show that there are larger biases in constraining SFRs and 
mean ages of galaxies using the SED-fitting method. As can be seen in these figures, the 
mean SFRs derived from the SED-fitting method systematically underestimate the intrinsic 
mean values by ~ 65 % for U-dropouts, ~ 58 % for B-dropouts and ~ 62 % for V-dropouts, 
while the stellar population mean ages are overestimated by about factors of two for all three 
sets of dropouts. 

The bottom rows of figure^ El and[7]show the results when the redshift is allowed to float 
as an additional free parameter in the fit - which is analogous to the case of color-selected 
LBG samples without spectroscopic (or pre-calculated photometric) redshift information. 

The bottom rows of the first columns of figure O [6l and [7| show that when we allow 
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redshift to vary as an additional free parameter, the differences between the distributions of 
SAM intrinsic and SED-derived stellar masses increase especially for B- and V-dropouts. The 
mean values of the SED-derived stellar masses are underestimated by 25 % (U-dropouts), 
51 % (B-dropouts), and 43 % (V-dropouts) when redshifts were allowed to float. 

More interestingly, bimodalities in the distributions of SED-derived SFRs and mean 
ages become more significant, as can be seen in the second and third columns of the bottom 
row of figures El and [71 The change is more significant for B- and V-dropouts. When we 
fix the redshift, we can only see hints of the existence of this bimodality in the SFR and 
age distributions of B-dropouts. The existence of these bimodalities, which are absent in 
the intrinsic distributions, indicates that there are sub-populations of LBGs whose behavior 
in the SED-fitting procedure is distinct from others, and the different behaviors of these 
subpopulations are exaggerated when we do not fix the redshifts. The directions of these 
bimodal distributions show that for this sub-population of LBGs, the SED-fitting method 
does not underestimate (or even overestimates) SFRs and underestimates mean ages. The 
characteristic SFHs of this sub-population of galaxies are discussed in § 14.41 

SED-fitting generally underestimates redshifts slightly and the ranges of redshift dis- 
crepancies, {zsED — zsam)/{^+zsam) are ~ -0.137 — 0.043 for U-dropouts, ~ -0.115 — 0.044 
for B-dropouts, and ~ -0.052 — 0.010 for V-dropouts. For galaxies whose redshifts are 
severely underestimated, stellar masses and ages are severely underestimated and SFRs are 
overestimated. The means of \zsed — zsams\/ i'^+^sAMs) are 0.022, 0.032, and 0.014 for U-, 
B-, and V-dropouts. 

Table. [2] summarizes the main results of the SED-fitting in the case where we fix the 
redshifts to the values in the SAM model galaxy catalog (analogous to observations with 
spectroscopic redshifts) as well as in the case where we allow redshift to fioat as an additional 
free parameter (analogous to observations without spectroscopic redshifts). The contents of 
table. [2] are the mean values of SAM intrinsic and SED-derived stellar masses, SFRs, and 
mean ages as well as redshifts for each set of dropouts. 

4.2. Biases in Constraining Physical Parameters 

In section 14. 1^ we showed the statistical behavior of various physical parameters derived 
from the SED-fitting compared with the intrinsic distributions. In this and the next section, 
we investigate more thoroughly the biases in the statistical properties of physical parameters 
derived from SED-fitting, focusing in particular on the dependencies of the offsets on various 
galaxy parameters. 
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First, in this section, we show how the discrepancies in LBGs' stellar mass, SFR, and 
mean age depend on the magnitudes and colors of LBGs. Specifically, we investigate the 
behavior of biases as functions of rest-frame UV and optical magnitudes (i.e. ACS Vqqq and 
IRAC "^4.5/im for U-dropouts, ACS i-rj^ and IRAC ms.g/im for B-dropouts, and ACS 2:850 and 
IRAC ^5.8/im for V-dropouts), and rest-frame UV and UV-optical colors (-B435 — Vme and 
^775 - "^3.6/im for U-dropouts, Veoe - ^775 and ^775 - ms.e/xm for B-dropouts, and ^775 - 2:850 
and 2850 ~ ''^4.5/im for V-dropouts) in figures [HlfTOl for the case of fixed redshifts and in 
figures [TTlfT^ for the case of allowing redshift as an additional free parameter. In these 
figures, relative errors of stellar mass, SFR, and mean age are defined as A^M* = {M^^sed — 
,sam)/M^:^sam, A'^SFR — {SFRsED — SFRsam)/SFRsam, and IS!' Age — {AgesED — 
AgesAM)/AgesAM, respectively. Here, valuesEO is stellar mass, SFR, or mean age derived 
from SED-fitting, and values am is the intrinsic stellar mass, SFR, or mean age of each 
galaxy. These figures clearly show that the SFRs of most galaxies are underestimated and 
mean ages are almost always overestimated. 

For both the case of fixed redshifts and where we vary redshift as a free parameter, the 
magnitude- and color-dependent behaviors are similar; however, the plots for B-dropouts 
when redshifts are fixed (figure [HD and for U-, B-, and V-dropouts when redshifts are allowed 
to vary freely (figures [mfT3|) reveal subpopulation(s) of LBGs whose behaviors in SED-fitting 
are distinguished from the majority of LBGs in each dropout sample even in stellar mass 
estimation. The bimodality which is clearly seen in figure [9] is not evident in the mass 
distributions shown in section 14.11 According to figure M and figures [TTlfTSl galaxies which 
are belong to this subpopulation show distinguished pattern of biases from the majority of 
galaxies: (1) Their stellar masses are more severely underestimated than the majority of 
galaxies. (2) SFRs are less severely underestimated or even overestimated. (3) Mean ages 
are underestimated for these galaxies in the subpopulation while they are overestimated for 
the majority of galaxies. Also, this subpopulation of galaxies is more likely to reside on the 
fainter side of the rest-frame UV magnitude distribution. A more detailed investigation of 
this subpopulation is given in section It should also be noted that figure [T2l indicates that 
there are actually (at least) two subpopulations whose behavior in the SED-fitting procedure 
are distinct from the majority of LBGs. 

The offsets in stellar mass, SFR, and mean age show relatively clear dependencies on 
rest-frame UV-optical color. (1) The stellar mass is more likely underestimated for redder 
LBGs, (2) the bluer the LBG is in rest-frame UV-optical color, the more the SFR is un- 
derestimated and the more the mean age is overestimated. Rest-frame UV-optical color is 
considered to be a crude indicator of a galaxy's stellar population age, so dependencies of 
offsets on rest-frame UV-optical color may indicate that (one of) the main cause(s) of the 
biases in SED-fitting is the mean age (or SFH), in the sense that systematic offsets between 
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intrinsic and best-fit SFRs and mean ages increase for younger galaxies. 

Interestingly, stellar masses tend to be overestimated for galaxies that are blue in rest- 
frame UV-optical color. These are also the galaxies whose biases in SFR and mean age are 
largest, which suggests that the mass-overestimation and age-overestimation share the same 
origin. 

The parameter that shows the clearest dependence on rest-frame optical magnitude is 
stellar mass. Stellar mass is more likely to be underestimated for LBGs which are brighter 
in rest-frame optical magnitudes, while stellar masses of fainter LBGs are more likely to 
be overestimated. At both faint and bright rest-frame optical magnitude, SFRs are clearly 
underestimated, while relative errors are distributed widely. This indicates that several 
different factors may contribute to the biases in SFR estimates. 

The offsets show almost no dependence on rest-frame UV-color, and relatively weak 
dependence on rest-frame UV magnitude. The dispersions of offsets are larger for galaxies 
whose rest-frame UV magnitudes are fainter, and SFRs are more likely underestimated for 
galaxies with bright rest-frame UV magnitudes. 

These color /magnitude dependencies become more complex when we do not fix the 
redshifts. More specifically, for some galaxies with red rest-frame UV-optical color, stellar 
masses are greatly underestimated, SFRs are severely overestimated, and ages are severely 
underestimated. These trends are similar for the subpopulation of B-dropouts shown in 
figure M 

Also, there are also some hints of galaxies with blue rest-frame UV-optical color, whose 
stellar masses are underestimated, SFRs are roughly correct, and ages are underestimated. 

Table. [3] lists the means and standard deviations of relative offsets for physical param- 
eters such as stellar masses, SFRs, mean ages, and redshifts for each set of dropouts. Here, 
the relative offset for each parameter is defined as (ValuesED — ValuesAAi) /{Value s am)- 
For redshift, the relative offset is defined slightly differently as {zsed — ^s'Aa/)/(1 + zsam), 
following the convention used in the literature. 

4.3. Origins of the Biases of Galaxy Population Parameters 

In this section, we investigate the dependencies of the fitting discrepancies on intrinsic 
properties of the SAM model galaxies, such as stellar mass, SFR, age, and specific SFR 
(SSFR; defined as SFR per unit stellar mass) to investigate the causes of these biases in 
SED fitting. 
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FiguresfT ^ [T5 l and [T6] show how relative errors (as defined in previous section) in stellar 
mass, SFR, and mean age are correlated with intrinsic stellar mass, SFR, mean age and 
SSFR for U-, B-, and V-dropouts, with redshift fixed. These correlations shed light on the 
origins of the biases found in sections 14.11 and 14. 2[ 

4-3.1. Origin of Bias in Age Estimation 

The clearest correlations are seen in 'A''"M=k vs. M^^sam\ 'A'^Age vs. AgesAM\ and 
IS!" Age vs. SSFRsam" ■ The tight correlation of relative mean age errors with intrinsic ages 
is not unexpected from the correlation between relative mean age errors and rest-frame UV- 
optical colors. This correlation and the one between relative mean age errors and intrinsic 
SSFRs are a strong indication that mean ages and/or SSFRs are the main cause of bias. The 
sense of bias is that the stellar population mean age overestimates are worse for galaxies with 
the youngest intrinsic ages and/or largest SSFRs (i.e. galaxies whose current SF activity 
is strong compared with the past SF activity), as can be seen in figure [T71 In this figure, 
which shows the intrinsic mean ages (y-axis) and SSFRs (x-axis) of B-dropout galaxies, blue 
dots represent galaxies whose relative age errors are 0.0 < IV Age < 0.75, i.e. galaxies 
with the smallest A'^Age. Green dots are for galaxies with 0.75 < A^Age < 2.0, and large 
red dots are for galaxies with the largest A^Age (> 2.0). This figure clearly shows that 
galaxies with very large age overestimates are relatively young galaxies with high SSFRs. 
The Young mean ages and large SSFRs of these galaxies indicate that they have experienced 
a relatively high level of SF activity recently. The mismatch between the SFHs predicted 
by the semi-analytic galaxy formation model and the simple SED templates from the BC03 
stellar population model is expected to be largest for galaxies with this type of SFH. In the 
BC03 templates, SFRs are assumed to decrease exponentially, so the strongest SF activity 
always occurs at early times, which is nearly the opposite of the SFHs of these galaxies. 
This difference makes ages overestimated severely for these galaxies. Figures [18] and [191 
which show SF activity as a function of lookback time for individual galaxies, support this 
speculation. Figure [TS] shows the typical SFHs (we will refer this type of SFH as 'type-2' 
from now on) of galaxies whose age-overestimation is largest. The star formation histories 
of these galaxies are clearly distinct from the SFHs of galaxies whose age-overestimation is 
smallest. The typical SFHs of galaxies with the smallest age-overestimation are shown in 
figure \W\ ('type-l' SFHs from now on). In figure [T71 there is a dearth of galaxies with old 
age and high SSFR and galaxies with young age and small SSFR. Very young galaxies with 
low SSFR would not pass our color selection criteria and/or magnitude limits or would not 
be detected in real galaxy samples. Old galaxies with high SSFR, in contrast, would be 
probably detected. To have a high SSFR and an old mass-weighted stellar population age. 
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galaxies would need to have had a very extreme SFH - for example, two strong, very short, 
widely separated bursts. 

4-3.2. Origins of Bias in Stellar Mass Estimation 

How does this age overestimation affect other derived physical parameters, such as stellar 
mass and SFR? When the mean stellar age is overestimated, some portion of the galaxy's 
luminosity will be attributed to older stars, with consequently higher mass-to-light ratios, 
than would be the case for the true SFH. This leads to an overestimation of the stellar mass. 

In previous studies, it has been suspected that the stellar masses of galaxies are under- 
estimated through single-component SED-fitting, since light from the recent star formation 
can easily mask some portion of the older stellar population. This is confirmed through the 
experiments done with the composite BC03 templates in section 15.41 As stated in section 
14. H the mean values of stellar masses are underestimated by ~ 19-25 %, which is in qual- 
itative agreement with the arguments of previous studies. However, figures [HlfTOl and [TUfTGl 
show that the stellar mass is not always underestimated. For some galaxies with very small 
stellar masses and/or relatively young ages, the stellar mass can be overestimated. These 
galaxies are the ones for which the age overestimation is large. We speculate that the stellar 
mass estimation of LBGs is affected by two different factors. One factor is the fact that 
the recent generation of star formation can dominate the broadband SED, leading to the 
underestimation of stellar mass. The other is the SFH difference between SAM galaxies and 
BC03 templates. If the age is overestimated due to the SFH difference, the mass-to-light 
ratio is overestimated, which results in the stellar-mass overestimation. 

In figures [TUfTGl it can be seen that the stellar masses are most likely underestimated for 
the galaxies with the oldest ages and/or smallest SSFRs, for which the age-underestimation 
is minimal. More clearly, (a)-(c) in figure [20] show that the stellar masses are underestimated 
when the age estimation is nearly correct. This confirms that the stellar mass is underes- 
timated in single- component SED-fitting when the effect of the other origin of bias - SFH 
difference between the SAM galaxies and the BC03 templates - is minimal 

In figure [20] (a)-(c), we can also see that the age discrepancies and the stellar-mass 
discrepancies are correlated. Stellar masses tend to be underestimated for galaxies with the 
smallest A'^Age, while they are overestimated for most galaxies with the largest A'^Age. 
This supports the speculation that the age-overestimation (due to the SFH mismatch be- 
tween the SAM galaxies and the BC03 templates) causes the mass-overestimation through 
mass-to-light ratio overestimation. Evidently, in the stellar mass estimation, two sources of 
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bias are compensating with each other. While SFH differences between the SAM galaxies 
and the BC03 templates tend to make stellar masses to be overestimated in SED-fitting, 
recent SF activity can easily mask older stellar population causing the stellar masses to be 
underestimated. The compensating effects of these two biases in the estimation of stellar 
mass explains why the stellar mass distributions are recovered better than other parameters 
through the SED-fitting, and also provides a clue as to why the stellar mass h as turned out to 



be the most robust pararneter in earlier studies based on the SED-fitting (e.g. lPapovich et al. 



20011 : IShaplev et al.ll200ll . [20051). 



4-3.3. Origins of Bias in Star Formation Rate Estimation 

Let us now consider the sources of bias in the SFR estimates. First, when ages are 
overestimated, the SED fitting erroneously assigns some portion of the luminosity to older 
stars instead of stars that are just forming. Figures (d)-(f) of figure [201 show that the SFRs 
are most severely underestimated for galaxies with the largest A^Age. 

However, the SFR still tends to be underestimated even when the age estimates are 
nearly correct. Even some of the most extreme underestimates can be found for galaxies 
whose SED-derived age is correct to within a factor of two. This indicates that there is 
another source of bias in the SFR estimation. 

Interestingly, the correlation between the relative SFR error and the relative stellar mass 
error shown in (g)-(i) of figure [20] reveal distinct behaviors of the upper envelopes in this 
correlation between galaxies with the positive relative stellar mass errors and galaxies with 
the negative errors, providing another indication that there are two sources of bias in the 
SFR estimation. For galaxies whose stellar masses are underestimated - i.e. galaxies whose 
age-overestimation is small due to the little SFH difference or galaxies with type-1 SFH - 
the SFR discrepancy is proportional to the stellar mass discrepancy. As stated in § 13. 2[ the 
SFR is calculated from the estimated stellar mass accumulated over 100 Myr. Thus, for 
galaxies with the smallest age discrepancy (with the type-1 SFH), the underestimation of 
stellar mass results in the SFR underestimation. 

However, unlike in the stellar mass estimation, both of these two origins of bias - 
the age-overestimation due to the SFH mismatch and the hidden old stellar population by 
the recent SF activity - cause the SFRs to be underestimated. This results in the overall 
underestimation of the SFR distributions and the large offsets in the mean SFRs. 

In addition to these two origins, the well-known 'age-extinction degeneracy' leads to 
more significant SFR underestimation. For galaxies with large (positive) A^Age, dust ex- 
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tinctions tend to be underestimated due to the 'age-extinction degeneracy'. This further 
deflates estimated SFRs for these galaxies. 

In summary, the main origins of biases in estimating physical parameters, such as stellar 
mass, SFR, and mean age, are: (1) the differences in the assumed SFHs in the SAM galaxies 
and in the BC03 stellar population templates, (2) the effects of the recent SF activity hiding 
some portion of old generations of stellar population, and (3) the age-extinction degeneracy. 
In the stellar mass estimation, issues (1) and (2) compete with each other, resulting in the 
best-fit stellar mass distributions that resemble the intrinsic distributions. For the SFRs, all 
of these issues work in the same direction, leading to the large offsets in the distributions 
and in the mean values. The mass- weighted mean ages are mostly affected by issue (1). 

Figures and figure [H] are similar plots with figures [TUfTBI and [201 but when red- 

shift is an additional free parameter in the SED-fitting procedure. These figures show trends 
similar to the redshift-fixed case except that subpopulations are more evident, especially for 
B- and V-dropouts. The increased ambiguity due to the lack of redshift information evi- 
dently enlarges the subsets of galaxies that behave distinctly from the majority of galaxies in 
the SED-fitting. However, for the majority of galaxies, the lack of redshift information does 
not significantly affect the SED-fitting results, which is not surprising given the relatively 
small mean redshift errors in SED-fitting. 

There is a subset of galaxies whose stellar masses are severely underestimated, SFRs 
are severely overestimated, and ages are severely underestimated. These are similar trends 
with those shown by a subpopulation of B-dropouts when redshift is fixed. However, the 
number of galaxies which belong to this subpopulation substantially increases when the 
redshift is allowed to float as an additional free parameter. For these galaxies, redshift is 
underestimated in the SED fitting. For U-dropouts, this sub-population is not as significant 
as for B- or V-dropouts. 

Also, there are galaxies with high SSFRs/young ages which act differently if we perform 
the SED-fitting without fixing redshift. For these galaxies, ages are greatly overestimated, 
stellar masses are overestimated, and SFRs are greatly underestimated when we fix redshift 
in the SED-fitting (i.e. these are the galaxies with type-2 SFHs). When we vary redshift 
freely as a free parameter, redshifts are slightly underestimated and ages/stellar masses are 
underestimated. The SFRs are similar to or slightly higher than the intrinsic values. 
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4.4. Characteristics of Sub-populations in the Fitted Distributions 

Figures [9l [T5| and (b), (e), and (h) of figure [201 reveal the presence of a subpopulation 
of B-dropout galaxies whose behavior in the SED-fitting is distinct from the majority of 
galaxies. For these galaxies, the ages are underestimated, the SFRs are overestimated, and 
the stellar masses are more severely underestimated than other galaxies. What makes the 
behavior of the galaxies in this subpopulation different from the majority of galaxies? 

Because one of the main origins for biases in SED fitting is the SFH difference between 
the SAM galaxies and the BC03 templates, it is plausible that the SFHs of these galaxies are 
distinct from others. Figure [25] shows SFHs of typical model galaxies in this subpopulation 
('tj^e-3' SFHs from now on). Generally, they have small SSFR values like those shown in 
figure [191 However, the SF activity in figure [251 shows a slower increase and more rapid 
decrease with a peak at later time compared with the galaxies shown in figure [191 The 
gradual decreases of the SFRs shown in figure [191 are not significantly different from the 
exponentially decreasing SFRs assumed in the BC03 model, which makes the relative age 
errors small for these galaxies. On the other hand, for galaxies shown in figure [251 the 
strong SF activity, which occurred relatively recently, dominates SEDs. Combined with 
the age-extinction degeneracy, this causes the mean ages to be severely underestimated, 
distinguishing behaviors of these galaxies in the SED-fitting. The purple crosses in figure 
[T71 represent galaxies in this sub-population. They are not clearly distinguished from other 
galaxies with the type-1 SFH in this age-SSFR domain, but have, on average, slightly younger 
ages than the type-1 SFH galaxies with similar SSFR. 

Of course, the SFHs of all galaxies are not clearly divided into typical examples shown 
in figures [TSl [THl or [231 For example, when only ACS and IRAC fiuxes are used (i.e. if 
smaller number of passbands are used; see § l5.2.2p . more galaxies behave similarly to galaxies 
with SFH type-3, making the bimodal distributions of SFRs and ages more prominent. 

5. Effects of Parameter Changes on the Results of the SED-fitting 

In the following sections, we investigate how the bias in the SED fitting behaves as we 
change some of the conditions in the SED-fitting procedure, such as the range of e-folding 
time of star formation history, r, combinations of broad passbands used, or the assumed 
SFH. 
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5.1. Two-Component Fitting 



Here, we try to allow more complex star formation histories in the SED-fitting, by 
using the two-component templates instead of single component ones. Some studies tried 
this method to constrain the hidden mass in old stellar population or for better estima- 
tion of total stellar mass. To construct the two-component stellar population templates, 
various combinations of simple SFHs have been tried in the literature, including: (1) com- 
bining a maxirnally old, instantaneous bu rst with a more moderately decreasing SFH (e.g. 
Papovich et al.l l200ll : IShapley et al.l 120051 ) . or (2 ) adding a secondary young bursty SFH 
component to an old , slowly decreasing SFH (e.g. iKauffmann et al.ll2003l : iDrory et al.ll2005l : 
Pozzetti et al.l 120071 ). Also, different authors have used different ways of fitting two com- 
ponents: (1) fitting with a young SF c omponent first, then fitting the residual SED with 



an old component (IShapley et al.ll2005l ). or (2) constructing the co mbined templates with 



various ratios betwe en a young and an old simple SFH templates ( IKauffmann et al.l 12003 
Pozzetti et aPboOTh . 



5.1.1. Slowly Decreasing Star Formation Histories with a Secondary Burst 

Here, we constructed the two-component SFH templates by adding (maximally) old, 
very slowly varying SFH templates (r = 15 Gyr), and younger, more bursty templates (r 
= 0.2 Gyr). Old components are assumed to start forming at Zf = 10. The star formation 
activity of the secondary burst is constrained to be initiated at least 200 Myr later than 
that of the old component, and at most 500 Myr earlier than the observed time, with the 
percentage of young templates varying between 5 % to 95 % (with 5 % step size). This 
construction is expected to reflect better the SFHs of some galaxies in the semi-analytic 
models used in this study (for example, galaxies with type-2 SFHs shown in figure [T5]l . 
and using this type of composite templates ought to reduce the systematic bias for galaxies 
with type-2 SFHs. The SFHs constructed in this way can mimic the SFHs of galaxies that 
experienced the secondary star formation due to merger/interaction. 

The trends in the SED-derived distributions of various parameters in this type of two- 
component fitting are: (1) the stellar masses are underestimated more severely than in 
the single-component fitting (middle row of figure [26]) . (2) the SFR/age distributions show 
reduced offsets compared with the single- component fitting (middle rows in figures [27| and 
l28l) . and (3) the bimodalities that existed in the age distributions have disappeared, while the 
bimodalities in the SFR distributions are enhanced for B- and V-dropouts (for U-dropouts, 
the SFR distribution becomes much broader). These bimodalities in the SFR distributions 
are, however, not driven by the subpopulation of galaxies with type-3 SFHs, as can be seen 
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below. 

The first two trends are not unexpected since we are adding a young stellar population 
with a low mass-to-light ratio to the SED templates. This young component makes the 
best-fit ages younger, thereby reducing bias in the mean ages, while its lower mass-to-light 
ratio makes the best-fit stellar masses smaller, thereby increasing bias. Younger ages, and 
thus increased dust extinctions, lead to higher SFRs, and decrease bias in the SFRs. 

Figure [21] - showing the correlation between the relative errors arising in the two- 
component fitting and the relative errors in the single-component fitting - reveals behaviors 
of galaxies with different SFH types in the two-component fitting performed in this section. 
Galaxies whose ages are greatly overestimated in the single-component fitting (with type- 
2 SFHs) show significantly reduced age-overestimation in the two-component fitting. This 
improvement in the age estimation is expected. In the single-component fitting, the ages are 
greatly overestimated for the SFH type-2 galaxies because the young component is ignored. 
By adding a young component in the templates, ages are better fitted for this type of 
galaxies improving the age estimation. For these galaxies, the stellar masses are generally 
overestimated in the single-component fitting. Improving the age estimates also improves 
the stellar-mass estimates 

However, when the age is relatively well constrained in the single-component fitting 
(i.e. for galaxies with type-1 SFHs), an added young component causes the underestimation 
of ages, and therefore more severe underestimation of the stellar mass which is already 
underestimated in the single-component fitting. The age underestimation of these galaxies 
couples with the age-extinction degeneracy, leading to the overestimation of the SFRs shown 
in (g), (h), and (i) in figure [291 This causes the bimodalities in the SFR distributions shown 
in middle row of figure [27[ 

Interestingly, for a small number of galaxies, the ages and the stellar masses derived in 
the two-component fitting are higher than the ones derived in the single-component fitting. 
These galaxies are the ones with SFH type-3 and are manifested as a subpopulation in 
figures M and [T5J Through the single-component fitting, ages and stellar masses of these 
galaxies are greatly underestimated, and the SFRs are severely overestimated due to the 
large age underestimation. The higher values of the ages and stellar masses derived through 
the two-component fitting reduce the errors in the age and stellar mass for these galaxies. 
For these galaxies with the type-3 SFHs, the best-fit ts are small in the single-component 
fitting. So, the actual effect of the two-component fitting performed in this section is to 
add an old component for these galaxies, while for the other galaxies (with the type-1 or 
type-2 SFHs), the the effect of the two-component fitting is to add a young component. This 
makes the best-fit ages from the two-component fitting older than the ones derived from the 
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single-component fitting, the best-fit masses higher (due to higher mass-to-light ratios), and 
the best-fit SFRs lower for galaxies with type-3 SFHs. Older best-fit ages of the type-3 
SFH galaxies derived in the two-component fitting removes the bimodalities shown in the 
SED-derived age distributions in the single-component SED-fitting. 

The age distributions derived in the two-component fitting are narrower than the in- 
trinsic distributions (middle row of figure [28]). The ages tend to be underestimated for older 
galaxies (possibly with SFHs of type-1 or type-2), while ages are more likely overestimated 
for younger galaxies (probably SFH type-3 galaxies). The SFR distributions derived in the 
two-component fitting are more extended toward the high SFRs than in the single-component 
fitting for U-dropouts (bringing the SFR distributions closer to the intrinsic distributions). 
For B- and V-dropouts, the bimodality in the SFR distribution is enhanced, i.e. the SFRs 
are overestimated for more galaxies. 

In summary, the two-component fitting performed in this section reduces bias in the SFR 
and age distributions, but increases the offsets in the stellar mass distributions. However, 
the detailed investigation reveals that the changes of behavior in the two-component fitting 
compared with the case of the single-component fitting are different for galaxies with different 
types of SFH. Errors in the estimation of ages and stellar masses are reduced for galaxies 
with SFHs type-2 or type-3, while the stellar mass errors increase for galaxies with type-1 
SFHs. 



5.1.2. Maximally Old Burst Combined with Slowly Varying Younger Components 

In the previous section, we experimented with two-component fitting by adding a young, 
burst-like (r = 0.2 Gyr) component to a more continuously varying (r = 15 Gyr) old 
component. Such two-component templates are expected to match better the type-2 SFHs, 
and turned out to give better age estimates for the galaxies with type-2 as well as type-3 
SFHs. 

In this section, we perform another type of two-component fitting in an attempt to 
give a better constraint on the hidden old stellar mass. To achieve this we construct the 
two-component SFH templates in a similar way done as in § 15.1.11 but exchange the roles 
of a r = 0.2 Gyr component and a r = 15 Gyr component. We add an old, r = 0.2 Gyr 
component formed at Zf = 10 and younger r = 15 Gyr components with various ages. The 
star formation activity of the young components is constrained to start at least 200 Myr 
later than an old burst, to make the two-component templates clearly distinguished from 
the single-component template. By adding an old, bursty component to a more continuously 
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varying SFH component, we can expect that this method will give us higher mass than the 
single-component fitting. 

The bottom row of figure [22] shows that the stellar mass distributions are moved toward 
higher values than the ones derived in § 15.1.11 (middle row of figure [26|l and also than the ones 
from the single- component fitting. Compared with the intrinsic distributions, the derived 
stellar mass distributions from the two-component fitting with an additional old, burst-like 
component are slightly crowded at high stellar mass for U- and B-dropouts. The bimodalities 
both in the SFRs and ages have disappeared (bottom rows of figures [271 and 1281) . 

The mean values of stellar masses are higher by ~ 19 %, 54 %, and 1 % than the values 
from the single-component fitting for U-, B-, and V-dropouts, respectively. For individual 
galaxies, the stellar mass from the two-component fitting with an old, burst-like component 
can be as large as several times of the stellar mass from the single-component fitting. For 
a few galaxies (mostly with type-3 SFHs), the stellar mass from the two-component fitting 
with an old burst can reach ~ 4 — 9 times of the stellar masses from the single-component 
fitting. However, the relatively small increase in the mean values of the best-fit stellar mass 
(especially for V-dropouts) indicates that the young component dominates the SEDs even 
in the two-component fitting performed in this section. 

The dominance of the young component can be seen by the fact that the mean values of 
the best-fit ages derived in the two-component fitting are much younger (by ~ 39 %, 35 %, 
and 42 % for U-, B-, and V-dropouts) than the mean values derived in the single-component 
fitting. (The star-formation time scale of the young component in the two-component fitting 
performed in this section is fixed as r = 15 Gyr, and this leads to younger best-fit ages - see 
§ [5X21 below.) 

Even though the two-component fitting performed in this section gives higher values 
of the mean stellar mass, it is not always true that stellar masses derived from the single- 
component fitting and this kind of two-component fitting bracket the true, intrinsic stellar 
mass. For some galaxies with small stellar mass (with log (M^/Mq) < 9.5), even the single- 
component fitting overestimates the stellar mass. On the other hand, the stellar mass derived 
through the two-component fitting with an old, burst-like component often remain smaller 
than the intrinsic value for some massive galaxies. 

If we were to limit the fractional contribu tion of the young compo nent to smaller values 



than allowed here - as done, for example, in [Kauffmann et al.[ (|2003[ ) - the stellar masses 



derived in the two-component fitting would become higher than the ones derived in this 
section. Also, we can derive higher stellar masses from the two-component fitting: (1) by 
fitting the old component first, then fitting the younger component to the residual fluxes (i.e. 
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forcing the contribution from the old component to increase), or (2) by setting the formation 
redshift {zf) of the old component higher (i.e. increasing the mass-to-light ratio of the old 
stellar component). However, even with these more extreme settings, it is still possible that 
the derived stellar masses for some very massive galaxies will be smaller than the intrinsic 
ones. 

In summary, it is not universally true that the single- component fitting and the two- 
component fitting (with an old, burst component added on more continuous SFH compo- 
nents) bracket the true stellar mass. 



5.2. Effects of Wavelength Coverage 

The main results presented in § 4 are based on the analysis using broadband photometric 
information from observed- frame optical through MIR range - i.e. ACS -B435- to Zgso-bands, 
ISAAC J- to Ks-hands, & IRAC 3.6 through 8.0 fim. However, not all the observed 
LBGs have photometric data with this wavelength coverage. Before the Spitzer era, the 
majority of the observed photometric data only covers up to the observed-frame NIR range. 
Thus, it is interesting and important to examine how the results vary as we use different 
combinations of passbands in the SED-fitting. 



5.2.1. SED-fitting without IRA C Data 



Several authors investigated the effects ( 
coverage of ~ 3-10 iJ,m) in constraining the i 


3f inclusion of IRAC photometry (0 
)roperties of high-redshift galaxies ( 


" wavelength 
Labbe et al. 


2005; 


Shaplev et al. 


2005; 


Wuvts et al. 


2007 


Eisner et al. 


2008 


)• 


Shaplev et al. (12005) anal- 



ysed z ^ 2 star-forming galaxies with and without IRAC photometric data. They reported 
that the SED-derived stellar mass distribution shows little change with the inclusion of IRAC 
data, while including IRAC data can reduce errors in the stellar mass estimation for indi- 
vidual galaxy. Investigating 13 2; ~ 2 — 3, red ( Js — Kg > 2.3) galaxies, iLabbe et al.l (120051 ) 
showed that the best-fit ages are younger without I RAC data for dusty star forming galaxies, 
while there is little change for old, dead galaxies. IWuyts et al.l (120071 ) studied 2 < z < 3.5, 
K-selected galaxies, and showed that inclusion of IRAC data does not chang e the overall 
distrib utions of stellar masses and ages. Analyzing < z < 5 observed galaxies, lElsner et al. 
( 120081 ) showed that the mean stellar masses increase when derived omitting Spitzer data. 
The discrepancy is maximum at z ~ 3.5 with log{Mu _i/ Mjj _k) ~ —0.5, and decreases with 
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redshift at 2 > 3.5. (At z < 3.0 or z > 4.0, log{Mu-i/ Mu-k) < -0.3. )0 In addition, they 
reported that despite this overall trend, the stellar masses and mean ages decrease without 
IRAC data for some very young (faint) galaxies. Thus, there is no clear consensus in the 
literature on the benefits of including IRAC. 

The effects of omitting IRAC data are shown in the third rows of figures [5l [6l and [7]). 
Without IRAC data, the stellar mass estimates shift to lower values for U- and B-dropouts, 
but are virtually unchanged for V-dropouts. Conversely, the SFR estimates are virtually 
unchanged for U- and B-dropouts, but shift to lower values (increasing the discrepancy with 
the SAMs) for V-dropouts. The age distributions shift to younger ages for the U- and B- 
dropouts (which actually brings them into better agreement with the intrinsic distributions) 
while for V-dropouts, the age distribution is only slightly changed. 

Evidently, removing the IRAC photometry for the U- and B-dropouts increases the 
dominance of the younger stellar populations due to the shorter wavelengths. This drives 
the best-fit ages to lower values, resulting in lower stellar masses for a given amount of 
stellar light. The stellar masses are reduced by ~ 57 % and ~ 48 % for U- and B-dropouts, 
respectively, relative to the results from fits, in which the IRAC data are included. 

The age and stellar mass decreases are largest for the galaxies with the high SSFRs 
and young ages, i.e. galaxies with type-2 SFHs (similar with the ones shown in figure [T8|) . 
These galaxies have roughly two components of stellar populations - an 'old', slowly varying 
component and a 'young', burst-like component. In the SED- fitting with IRAC photometry 
included, the 'old' component dominates the SED, resulting in much older best-fit ages than 
the intrinsic ages. However, without IRAC photometry, the SEDs cover only up to the rest- 
frame ~ 4000 — 5000 A. Due to the resulting shortage of information at long wavelengths, 
the 'young' component comes to dominate the SED. This moves the best-fit ages younger, 
even younger than the intrinsic ages in extreme cases. 

The changes of the best-fit age and stellar mass for V-dropouts are much smaller (~ 13 
% and ~ 4.8 % of decreases, respectively). This is presumably due to the generally younger 
ages of V-dropouts. As can be seen in figures dS] and dni V-dropouts, on average, started 
forming stars more recently than U- and B-dropouts. Therefore, the proportion of the old 
stellar population hidden without IRAC data is much smaller for V-dropouts than U- and 
B-dropouts. 

Interestingly, the bimodalities in the SFR- and age-distributions shown (for B-dropouts) 
in figures [3] and m disappear when we exclude IRAC photometry. For the type-3 SFHs shown 



'^Mu^i and Mjj-k are stellar masses that are derived with and without IRAC data, respectively. 
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in figure [251 the SFRs have lower values and the ages (and stellar masses) have higher values 
than the ones derived including IRAC bands. The behaviors of these galaxies are thus 
opposite to the majority of B-dropouts. 

What makes this sub-population of galaxies (with type-3 SFHs) behave differently from 
other galaxies? The difference of the SFHs between the ones shown in figure [25] (i.e. type-3 
SFHs) and the ones shown in figure [19] (i.e. type-1 SFHs) becomes significant when the 
lookback time is larger than ~ 400 — 500 Myr. With only ACS and ISAAC photometry, 
which covers only up to the rest-frame ~ 4000 — 5000 A, the SFHs at early time are hard to 
constrain. Therefore, the SED-fitting without IRAC data cannot discriminate between SFH 
type-1 (figure [T9]) and type-3 (figure [25]) . resulting in the disappearance of the bimodalities 
in the SFR and age distributions. 

Figure [3D] shows the ratios of the best-fit stellar masses, SFRs, and ages with and without 
IRAC photometry for U-, B-, and V-dropouts. Without IRAC photometry, the stellar masses 
are underestimated (compared with when IRAC data are included) for most U-dropouts 
(figure [5U]-(a)). and for the majority of B-dropouts (figure [5D]-(b)). For the subpopulation 
of galaxies in B-dropouts, the best-fit stellar masses without IRAC data are larger than the 
ones with IRAC data. For V-dropouts (figure [30]- (c)). the stellar masses are underestimated 
for some galaxies, and overestimated for other galaxies, making the distributions with and 
without IRAC data similar. 

In summary, the effect of removing the IRAC data depends on the redshift and/or 
SFHs. This redshift and SFH dependence can explain the apparent disagreement between the 
different previous investigations, since the samples included galaxies with different redshift 
ranges and also different types of galaxies with possibly different SFHs. 

5.2.2. SED- fitting without ISAAC data 

Next, we perform the SED-fitting with only ACS and IRAC photometry, excluding 
the J, H, and i^"s-band ISAAC photometry. The effect of omitting ISAAC photometry is 
insignificant for U-dropouts and for the majority of B- and V-dropouts. As can be seen by 
comparing the second and the fourth rows from the top in figure [5] the distributions of best- 
fit stellar masses, SFRs, and ages show little change without ISAAC data for U-dropouts. 
The mean stellar mass and mean age increase by 1.9% and 2.4% compared with the values 
derived using ACS + ISAAC + IRAC photometry. 

For B- and V-dropouts, the bimodalities in the SFR/age distributions become more 
prominent as can be seen in the fourth rows of figures [H] and [7] For some galaxies whose age 
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offsets are small when the SED-fitting is done with full photometry, the best-fit ages become 
much younger if we use only ACS and IRAC photometry. This age underestimation leads 
to the stellar mass underestimation and SFR overestimation as explained in § 14.3^ and the 
extinction overestimation due to the well-known age-extinction degeneracy enhances the SFR 
overestimation further. These galaxies skew the mean values of stellar mass and mean age 
of total sample lower by about 10 % for B- and V-dropouts, even though the majority of 
galaxies show little change. 

5.3. Effects of r-range Used in SED-fitting 

The changes of the allowed ranges of parameters, such as r, t, or metallicity, in the 
SED-fitting would affect the derived values of physical parameters, as well. Here, we focus 
on the effects of the different range of rs used in the SED-fitting on the estimation of physical 
parameters of LBGs. This investigation is beneficial for the comparison with previous works 
done using the SED-fitting methods with various r ranges as well as for better understanding 
biases of the SED-fitting methods. 

5.3.1. SED-fitting with t < 1.0 Gyr Templates 

First, we limit the r range to < 1.0 Gyr during the SED-fitting. Through this experi- 
ment, we can look into the biases arising due to the usage of 'not-long-enough' r values in 
the SED-fitting. 

As expected, the enforced smaller r values cause smaller best-fit t, to match t/r val- 
ues, compared with the case when the full range of r is allowed from 0.2 Gyr to 15 Gyr. 
This bias systematically makes the best-fit ages to be younger. As explained in § 14.31 this 
age- underestimation leads to the mass-underestimation, increasing differences between the 
intrinsic- and SED-derived stellar masses. The effects of the limited r values on the SFR 
estimation is complicated due to the age-extinction degeneracy. The offsets due to the lim- 
ited r values as < 1.0 Gyr are greatest for B-dropouts because of the severely enhanced 
bimodalities. Relative changes of the mean stellar masses and ages when we limit r range 
to < 1.0 Gyr, compared with the mean values derived with the full range of rs, are shown 
in table. IH 

As can be seen in figure EH the lowered best-fit values of ages/stellar masses are largely 
driven by galaxies with preferentially young intrinsic ages (top row) and/or high intrinsic 
SSFRs (middle row). These are mostly the galaxies with SFH type-2, whose ages are severely 
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overestimated when we fit the SEDs with the full range of rs (from 0.2 Gyr to 15 Gyr) 
(bottom row of figure [31]) . 

The type-2 star formation histories have two main components. One component is a 
relatively low level of long-lasting SF activity, which corresponds to large r values. The 
other is a relatively young, strong SF activity, which is more likely represented by small 
T. When we fit SEDs with full range of rs, best-fit models tend to be determined by the 
underlying, long-last SFH component giving the severely overestimated ages to these galaxies 
(see § 14.3.11) . In the case when only small values of r (< 1.0 Gyr) are allowed in the SED- 
fitting procedure, the best-fit models are more likely determined by the young SF component 
with relatively small star formation time scale. The best-fit t values then are much smaller 
than the ones derived utilizing the full range of rs, resulting in much younger best-fit mean 
ages and much smaller best-fit stellar masses. 

Another significant feature is the more prominent bimodalities in the distributions of 
best-fit SFRs and ages. Galaxies in the smaller sub-population have very young best-fit 
ages and very high best-fit SFRs compared with the remaining galaxies. Figure [HI] shows 
the increase in the number of galaxies which belong to this sub-population as an effect of 
limitation on the allowed r values. They are the galaxies with small intrinsic SSFRs (middle 
row) as galaxies with type-1 or type-3 SFHs. The discrepancies between the intrinsic- and 
best-fit ages are very small when the full range of rs is used in the SED-fitting (bottom row), 
which means they behave like the galaxies with type-1 SFH. However, when we restrict r 
as < 1.0 Gyr, their best-fit ages become much younger and join the sub-population. The 
bottom row of figure [31] also shows that galaxies which are in this subpopulation when r has 
the full range remain in the subpopulation when r is restricted to be < 1.0 Gyr. 

In summary, if we limit r to be < 1.0 Gyr in the SED-fitting, the overall trends are: (1) 
the best-fit ages (and hence the best-fit stellar masses) become smaller for the galaxies with 
type-2 SFHs. (2) a larger number of galaxies joins the subpopulation that has much smaller 
best-fit ages and stellar masses than the intrinsic SAMs values. These trends result in (1) 
an overall downward shift of the age/stellar mass distributions (increasing the discrepancies 
between the intrinsic- and SED-derived stellar mass distributions) and (2) more prominent 
bimodalities in the SFR/age distributions. 

5.3.2. SED-fitting with t=15 Gyr Templates 

Next, we hold r fixed at 15 Gyr, which is equivalent to assigning a constant star forma- 
tion rate, considering the age of the universe at redshifts ~ 3 — 5. Figure WI\ reveals that the 
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best-fit values of t^is generally increase for galaxies with small best-fit taii derived with full 
range of rs, but decrease if the best-fit taiis are large. Here, t^i^ refers the best-fit t derived 
if we fix r at 15 Gyr, and taii is the best-fit t obtained when we allow the full range of r. 
This behavior is caused by the restriction on t to be younger than the age of the universe 
at each redshift If the best-fit t is already large for the full range of rs, there is no room 
to increase t to match the red color of these old galaxies. Instead, the fitted value of the 
extinction increases and the fitted age generally decreases. 

Also, with the larger value of r, larger t does not always results in older mean ages while 
smaller t always makes mean age younger, since mean age is a function of r as well as t. 
Therefore, the mean ages (and the stellar masses also as a result) are slightly underestimated 
overall. Galaxies whose best-fit ts are very small (< 0.2 Gyr) show little differences in best-fit 
t with or without the r = 15 Gyr restriction. 

Table. H] shows the relative changes in the mean values of stellar masses and ages when 
we set r = 15 Gyr compared with the case when we allow the full range of r. 

5.4. SED Models with Extreme SFHs from BC03 Model 

We further test what would be the results of the SED-fitting for the galaxies with the 
extreme star formation histories (SFHs). To do this, here, we construct three types of toy 
models from the BC03 model, replacing SAM galaxies. The parameter settings used in these 
toy models are summarized in table. [51 The aims of each toy model are to examine the biases 
which arise: (1) when we use shorter rs than real in the SED fitting (toy model 1), (2) when 
we use longer rs than real (toy model 2), and (3) when we try the single-component fitting 
for the galaxies with clearly distinct, two generations of star formation. 

5.4-1 ■ Effects of SED-fitting Using Too Small ts 

First, in the case (1), the toy model SEDs have a very long SF time scale (r = 15.0 
Gyr). By restricting r for the SED fitting not to exceed 1.0 Gyr, we can re-examine more 
transparently (because we compare the same BC03 models) what would happen if we fit the 
SEDs with r values much shorter than the actual SF time scales of (model or real) galaxies. 

As expected, the best-fit mean ages are underestimated. To match colors (or SED 
shapes) of 'r = 15 Gyr' samples with much shorter rs, the SED templates with smaller te 
are found as the best-fit templates, which causes the mean ages are systematically underes- 
timated. The amount of the age underestimation increases with age. The relative age error 
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reaches up to ~ 30 % underestimation for the oldest galaxies (with t = 1.0 Gyr). 

The systematic underestimation of ages leads to systematic mass underestimation, since 
the given amount of light from galaxy is attributed to younger (more massive) stars with 
lower mass-to-light ratios. The relative mass underestimation increases as the relative age 
underestimation increases. For the oldest galaxies whose relative age underestimation reaches 
up to ~ 30 %, relative mass underestimations are ~ 9 - 19 %. The 'Younger-than- input' 
best-fit ages result in higher average SFRs, whose relative overestimation spans from ~ 
through 9 %. 

For the samples with t as small as 0.01 Gyr, all the physical parameters, such as stellar 
masses, SFRs, and mean ages, are well recovered through the SED-fitting. This reflects the 
fact that the effects of SFH difference are insignificant for the very young galaxies. 

In summary, if one tries to fit galaxies' SEDs with relatively short rs, the resultant 
best-fit stellar masses and mean ages can be underestimations of the true values especially 
for galaxies which have very extended SFHs for relatively long timescales. These trends 
confirm the speculation of § 15.3. 1[ The SFRs are overestimated, but the relative errors are 
not as large as those of stellar masses/mean ages. 

5.4-2. Effects of SED-fitting Using Too Large ts 

In the case (2), the input toy models have shorter r than the values allowed in the SED- 
fitting procedure. This highlights the biases that can arise in the SED-fitting if galaxies 
have much shorter SF time-scales (probably burst-like) than the r values allowed in the SED 
fitting. 

The direction of bias in the best-fit mean ages is divided into two regimes depending on 
the actual age (or t) of each galaxy. For the model galaxies with small enough t (i.e. t = 
0.1, 0.2 Gyr) compared with the age of the universe at corresponding redshift (which is 2; ~ 
4, in this case), the best-fit ts are larger than the input ts to match the SEDs with longer rs 
than the input (r = 0.2 Gyr), leading to the overestimation of the mean ages by amount of 
~ 30 % for the galaxies with the input t = 0.1 Gyr, and ~ 50 % for the ones with the input 
t = 0.2 Gyr. 

For the similar reason as in the case (1) (but, in the opposite direction), the age over- 
estimation results in the mass overestimation by attributing light to older stars with higher 
mass-to-light ratios. The mass overestimation is about ~ 11-13 %. The 'older-than- input' 
best-fit ages make the SFRs underestimated by ~ 22-34 %. 
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However, for the model galaxies with relatively large values of t, which are not much 
shorter than the age of the universe, there is not so much room for t to increase. So, for the 
toy model galaxies with input t = 1.0 and 1.3 Gyr, the best-fit ts are only slightly larger 
than the input ts, and the mean ages arc younger than the input due to larger rs than 
input. Here, the stellar masses are overestimated even for the underestimated mean ages, 
mainly because of the more extended SFHs. The important difference between the models 
with large input ts and small input ts is that the derived dust extinctions, parameterized as 
E{B — V), are greatly overestimated to compensate the 'younger-than- input' mean ages for 
the SED models with large input ts. 

During the SED-fitting procedure, the best-fit ts generally move in the direction which 
'correct' the difference in rs - 'smaller-than-input' ts for the long input r models, and 'larger- 
than-input' ts for the short input r models. However, for the toy model galaxies with small 
TS and large ts, the SED-fitting cannot overcome the r difference by adjusting the best-fit 
ts due to the restriction that t be smaller than the age of the universe. Instead, the best-fit 
E{B — V) starts to be overestimated by the amount of AE{B — V)^ 0.25-0.3 for the input-t 
= 1.0 Gyr models and ~ 0.45 for the input-t = 1.3 Gyr models. This large bias, in turn, 
results in a large bias toward overestimated SFRs. The best-fit SFRs arc about ~ 13-17 
times of the input SFRs for input-t = 1.0 Gyr models, and reach up to ~ 70-110 times of the 
input SFRs for the input-i = 1.3 Gyr models. The extinction overestimation also leads to 
the larger mass overestimation, but the effects are not as dramatic as in the SFR estimation. 
This example illustrates how the well-known 'age-extinction' degeneracy affects the results 
of the SED-fittings, especially for the galaxies with extreme SFHs and/or for the case when 
parameter (for example, r) space allowed during the SED-fitting is not sufficiently large. 

In summary, if one tries to fit galaxies' SEDs with very long rs only, the resultant stellar 
masses are generally overestimated. If galaxies have sufficiently young ages (compared with 
the age of the universe at the redshifts where galaxies reside) the mean ages are overesti- 
mated, as the best-fit ts to be much larger than the input ts to compensate the r difference. 
However, for old galaxies whose input ts are compatible to the age of the universe at cor- 
responding redshift, the direction of bias in the age estimation is that the mean ages are 
underestimated (even though the best-fit ts are slightly larger than the input ts). Instead, 
the dust extinctions are greatly overestimated, which makes the best-fit SFRs erroneously 
high (up to two orders of magnitude). 



-34- 



5.4-3. Effects of Single Component SED-fitting for Galaxies with Two Generations of Star 

Formation 

Lastly, in the case (3), the model galaxies have two clearly distinguished generations 
of star formation with t = 0.1 & 1.0 Gyr. Both of the components are set to have r = 0.2 
Gyr. When we try to derive the best- fit physical parameters of these two component model 
galaxies via the single-component SED-fitting, both the mean ages and stellar masses are 
underestimated, as the result of the older stellar population being at least partially ignored. 

The age underestimation is minimal (~ 10 %) for the toy model galaxies in which the 
young-to-old ratio is smallest (i.e. young component fraction ~ 0.1). The age underesti- 
mation increases as the proportion of young component increases (up to ~ 77-88 % for the 
ones with young /old ~ 1.0). This implies the best-fit parameters tend to be determined by 
the stellar component whose fractional occupation is large. However, as the young compo- 
nent fraction increases further to young/old = 2.6, the age underestimation decreases and 
becomes similar to that for the galaxies with young /old = 0.52, because the mean ages of 
these galaxies (with young /old = 2.6) are already small enough to be greatly underestimated. 

As stated in previous sections, the age underestimation propagates to the mass under- 
estimation. The mass underestimation is minimal (~ 9 — 11 %) for the model galaxies with 
young / old = 0.1 due to the lowest age underestimation, and for the ones with young /old = 
2.6 because they have the smallest portion of their mass in the old component. This means 
the mass underestimation depends on two factors - the degree of age underestimation (i.e. 
the amount of the mis-interpretation of the mass-to-light ratio) and the fraction of stellar 
mass in the old component. The mass underestimation for galaxies with young /old — 0.26, 
0.52, and 1.0, is in the range 21-40 %. 

The SFRs are underestimated by ~ 14-43 %, and the discrepancy is larger for the galax- 
ies with larger fractional mass in the old component. This is because the best-fit parameters 
are affected by the old component while the SFR contribution of the old component is very 
small (~ 0.2-4 % of the total SFR, depending on the mass fraction of the old component). 
The best-fit rs are larger than the input (i.e. 0.2 Gyr), and the best-fit ts are larger than 
the input young component (i.e. 0.1 Gyr) due to the effects of the old component, and this 
leads to lower SFRs. 

In summary, if one tries to fit galaxy's SED with single-component SED-fitting when the 
real SFH has two episodes, the mean ages and the stellar masses are generally underestimated 
indicating the old components are to some extent masked. And, the SFRs are underestimated 
since the old components with very little current SF activity pollute the SEDs. 



- 35 - 



6. SFR Estimation from Rest Frame Ultra-violet Luminosity 

For galaxies with a roughly constant star-formation rate, the extinction-corrected UV 
luminosity is expected to provide a reasonably good estimate of the star-formation rate. 
This is fortunate, because often the only data available for high-redshift galaxies are a few 
photometric data points in the rest-frame UV. 

The relation betwee n SFR and U V lum inosity can be calibrated using spectral synthesis 



models, such as BC03. iKennicuttI (119981 ) notes that the calibrations differ over a range 
of ~ 0.3 dex, when converted to a common reference wavelength and IMF, with most of 
the difference reflecting the use of different stellar libraries or different assumptions about 
the star- formation timescale. The calibrations usually assume constant or exponentially 
declining star-formation rates. 

It is interesting to see how well this technique works for the more varied star-formation 
histories of the semi-analytic models. In this case, we are using BC03 and the same IMF for 
both the calibration and the SAM galaxies, so the discrepancies in derived SFRs must be 
primarily due to the different star-formation histories. 

Here, we derive SFRs of SAM B-dropout galaxies with the assumption that we do not 
know their redshifts, which is similar to the case when there is no spectroscopic redshift 
information for a color-selected LBG sample. 

The (dust-uncorrected) rest-frame UV (Aq = 1500 A) luminosity of each SAM B-dropout 
galaxy is calculated from the ^775 band flux as 



^!^,1500 



+ 

assuming all B- drop outs are at z = 4.0. 



X fu,i775, (12) 



Here, fu,i775 and L,^ i^qq are specific flux at ^775 band and speciflc luminosity at rest-frame 
1500 A, respectively, and is luminosity distance at redshift z, which is assumed to be 4.0. 

There are two major possible sources of systematic biases, which can arise due to the 
assumption that all galaxies are at z = 4.0: (1) ignoring the bolometric correction - which 
is a decreasing function of reds hift and is small - a factor of 1.0 ± 0.2 in the redshift range of 



B-dropouts (3.6 < z < 4.6) for IChabrierl (120031 ) IMF and solar metallicity, and (2) ignoring 
error in luminosity distance (di). 

These two sources of bias act in opposite direction. Not including bolometric correction 
causes the UV luminosity to be slightly overestimated at high redshift, since it is a decreasing 
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function of redshift. Underestimation of the luminosity distances for galaxies at high redshift 
{z > 4.0) leads to underestimated UV luminosities. The correction factor due to the error 
in di estimation is slightly larger at high-redshift than the bolometric correction factor, 
and ranges from 0.8 to 1.4 in the redshift range of B-dropouts. Together, this will cause 
rest-frame UV luminosity to be slightly underestimated at the high-redshift end of the range. 

For all B-dropout galaxies, dust-extinction is assumed to be E{B — V) = 0.15, which is 



thesame val ue used for dust-correc t ion in high-redshift LBG studies, such as lGiavalisco et al. 
J2004bh and lSawicki fc Thompsor] J2006h . 



Star-formation rates are then calculated using the conversio n of iKennicuttI (Il998[) di- 



vided by 2.0 t o correct for the different assumed IMF, because iKennicutt 
Salpeterl JlQSsj ) IMF with a mass range 0.1 to 100 Mq. 



199 



uses a 



Figure [33] shows the ratio of dust-corrected, UV-derived SFR to intrinsic SFR as a 
function of redshift. In this figure, we can see the redshift dependent behavior of the star- 
formation rate calibration from the rest-frame UV. This behavior is a combined effect of 
luminosity-distance underestimation, which is small, though, as explained above, and the 
redshift-dependent difference in average dust-extinction. The large scatter of SFR ratios at 
a given redshift reflects the variation of galaxies' intrinsic dust-extinction and SFH. 

Differences in dust extinction can have a large effect. For example, small change in 
assumed value of mean E{B — V) will significantly change the derived SFR values - for a 
range of E{B — V) from 0.10 to 0.20, the dust-correction factor vary from 2.6 to 6.5. 

With the assumed dust-extinction of E{B — V) = 0.15, mean UV-derived SFR is 11.099 
Mq yr~^. For comparison, mean values of intrinsic and SED-derived SFRs are 15.650 
Mq yr~^ and 6.638 Mq yr~^, respectiveljlf]. 

Figure [31] shows the distributions of intrinsic SFRs (left), and of SFRs derived by 
assuming all galaxies are at z = 4.0 and assuming dust-extinction of E{B — V) = 0.15 
(right). The distribution of SFRs derived from rest-frame UV luminosity, assuming all 
galaxies are at z = 4.0, is much narrower than the intrinsic one. 



^ Since calibra tion between UV luminosity and SFR is derived assuming constant SFR over 100 Myr 
(|Kennicuttlll998i ). comparison between UV-derived SFR and intrinsic or SED-derived SFR, which is averaged 
over last 100 Myr, is relevant. 
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7. Discussion 

We have shown in this paper that there are significant biases in the physical parameters 
derived from the SED-fitting of standard r-model to broadband photometry. Biases are 
severe especially for the SFR and mean age, but even the SED-derived stellar masses are 
biased. We now address in some detail how the biases in the derivation of these physical 
parameters can affect the investigation of high- 2; galaxies and the inferred galaxy evolution 
studies. 



7.1. Artificial Age bimodality 



Figures [3] and S] show that there are bimodalities in the SFRs and mean ages derived 
through the SED-fitting (most clearly shown for B-dropout LBGs) while there is no such 
bimodality in the intrinsic distributions. These bimodalities are enhanced when we try to fit 
the SEDs omitting some available input information, such as NIR photometry from ISAAC or 
the spectroscopic redshift (figures [5l[6l and [7]). The main origin of these artificial bimodalities 
is the mismatch of the SFHs between the SAM model galaxies and the templates from the 
BC03 stellar population model as explained in § 14.41 Such bimodalities can lead to the false 
interpretation that there are clearly distinguished populations among the similarly selected 
star-forming galaxies. 



Recently, iFinkelstein et al.l (120081 ) analysed 14 Lyman-a emitting galaxies (LAEs) and 
found that there is a clear bimodality in their age distribution, which are derived through 
SED-fitting, such that their ages either very young (< 15 Myr) or old (> 450 Myr). Based 
on this bimodality, they concluded that there are two distinct populations of LAEs - dusty 
starbursts and evolved galaxies. However, acco rding to the resul t s pres ented in our work, it 
is possible that this age bimodahty reported in IFinkelstein et al.l (120081 ) may not be real but 
an artifact which arises in the SED-fitting procedure due to the difference of their SFHs. Of 
course, caution should be applied in interpreting the SED-fitting results of the LAEs based 
on our analysis performed for the LBGs. It is still controversial how similar (or how different) 
the LBGs and L AEs are, despite some indications of similarities in th eir physical parame ters 
( iLai et al.l 120071 ) and the possible overlap between LBGs and LAEs (jRhoads et al.ll2008l ). 



Shapley et al.l (120051 ) also reported a subset of galaxies in their sample of 2; ~ 2, star- 



forming galaxies with extremely young ages (with t < 10 Myr). The galaxies in this sub- 
population a re relatively le s s mas sive (<log M^/Mq> = 9.64 (calculated from column 7 
in table 3 of IShapley et al.l (120051 )). while the mean value of total sample is 10.32. More 
dramatic difference between galaxies in this subset and other remaining galaxies is shown 
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in the derived SFR distribution. Among 72 galaxies, there are 10 galaxies whose SFRs are 
larger than 200 Mq yr~^, and among these 10, nine galaxies belong to this subset of galaxies 
with t < 10 Myr according column 8 in their table 3 (one remaining galaxy with high SFR 
has the best-fit t = 15 Myr). In contrast, almost half of their sample has very low SFRs 
(< IOM0 yr^^). Based on our analysis, it is plausible that the 10 galaxies with the very 
high SFRs have similar SFHs as shown in figure [25] (i.e. type-3 SFHs), and therefore their 
ages and stellar masses are underestimated while their SFRs are greatly overestimated. If 
so, the ages and SFRs of their total sample would present more continuous distributions. 



Bimodalities in the inferred ages are also see n among z ^ 5 LBG s of IVerma et al.l (120071 ) 
and among 14, 2; ~ 5 LBGs with spectroscopy of IStark et al.l (120071 ). 



More interestingly, in the IShapley et al.l (120051 ) sample, there are three galaxies whose 
physical parameters derived through the SED-fitting methods do not agree with the indi- 
cations from their rest-frame UV spectra. The SED-fitting results indicate that these are 
old [t/r » 1) galaxies with very low level of SFRs (~ a few Mq yr~^), while there are 
indicatio ns of a young p opula tion of stars with active star formation in there spectra (see 
§ 4.5 in IShapley et al.l (120051 ) for detail). This apparent disagreement can be explained if 
these galaxies have the star formation histories similar with the ones shown in figure [18] (i.e. 
type-2 SFHs). 



7.2. Possible Descendants of High Redshift Star Forming Galaxies 



There is great interest in connecting LBGs to their possible descend ants - i.e. deter 



mining whether or not they are progenitors of local massive ellipticals (e.g. 



19971 : lAdelberger et al.lll998l : ISawicki fc Yeelll998l : ISteidel et al.lll998l : ISomerville et al 



Lowentha 



et al. 



200l|). 



Many properties of LBGs - including their sizes, morphologies, number densities, clustering 
properties, and physical properties - are closely linked to this issue. Therefore, more accu- 
rate estimation of LBGs' physical parameters are clearly important for addressing this issue 
properly. 

According to our analysis, the stellar masses derived using the single-component SED- 
fitting methods tend to underestimate the true values, on average, by 19-25 %. Moreover, 
bias in the stellar mass estimation seems to strongly depend on the stellar mass itself, in 
a sense that the stellar masses are more severely underestimated for more massive galaxies 
(see figures mi [151 and [T6l). For some of very massive galaxies (M^, > lO^^M©), the best-fit 
stellar masses can be less than half of the intrinsic stellar masses. Accompanied by the 
SFR underestimation, which is generally more severe than the stellar mass underestimation. 
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this discrepancy can significantly affect the discrimination among possible evolutionary de- 
scendants of massive LBGs. Since the current (or recent) SFR is a measure of the possible 
additional stellar mass which can be added to LBGs during their evolution to the lower 
redshift, the underestimation of the stellar mass and SFR of massive LBGs can plausibly 
lower their possible stellar masses at lower redshift or at 2; ~ significantly. 

Interestingly, studies based on the clustering properties of LBGs spe culated that mas 
sive, high-^: LBGs could be the p rogenitors of lo c al ma ssive ellipticals (e.g. lAdelberger et al. 
19981 : ISteidel et al.lll998l ). while ISawicki fc Yed (Il998l ). based on the SED-fitting analysis, 
suggested that LBGs would not become sufficiently massive at low-z to be massive ellipti- 
cals unless experience significant number of mergers. 

Biases and uncertainties in the estimation of LBGs' ages - which are shown not only to 
be large and but also to vary significantly as the fitting parameters change in our study - 
can also affect understanding of LBG properties and the predictions regarding their possible 
evolutionary paths. For example, biases in the age estimation would propagate into errors in 
their estimation of the star-formation duty cycle, which would, in turn, affect the estimates 
of the number of galaxies which have similar stellar masses/ages with detected LBGs, but 
are undetected due to their reduced SFRs. 



8. Summary & Conclusions 

In this paper, we examine how well the widely-used SED-fitting method can recover 
the intrinsic distributions of physical parameters - stellar mass, SFR, and mean age - of 
high-redshift, star-forming galaxies. To this end, we construct model high-redshift galaxies 
from the semi-analytic models of galaxy formation, make a photometric catalog via the BC03 
synthesis model, and select LBGs through the appropriate color selection criteria based on 
their broadband colors. Then, we perform SED-fitting analysis, comparing the photometric 
SEDs of these model galaxies with various galaxy spectral templates from the BC03 stellar 
population model to derive the distributions of best-fit stellar masses, SFRs, and mean ages. 
We use this test to explore (rather exhaustively) the errors and biases that arise in such SED 
fitting and the underlying causes of these errors and biases. 

Here are the summary of the main results of this work. 

1. When we fix the redshift to the given value in the SAM catalog and use ACS /ISAAC / IRAC 
passbands, the SED-fitting method reproduces relatively well the input distributions of stel- 
lar masses with a minor tendency to underestimate the stellar masses and with substantial 
scatter for individual galaxies. The mean stellar masses are underestimated by about 19~25 
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%, which is due to the fact that the old generations of stars can be hidden by the current 
generation of star formation. The distributions of SFRs and mean ages show larger offsets 
than the stellar mass distributions. The SFRs are systematically underestimated and the 
mean ages are systematically overestimated, and these trends mainly reflect the difference in 
the SFHs predicted by the semi-analytic models and assumed in the simple galaxy templates 
used in the SED fitting. The well-known 'age-extinction degeneracy' plays an important role 
in biasing the derived SFRs. 

2. When we use redshift as an additional free parameter, the discrepancy between the 
intrinsic- and SED-derived stellar mass distributions increase (i.e., the overall stellar mass 
underestimation becomes worse), while the bimodalities which appear in the SFR & mean 
age distributions become more significant. The distributions of offsets of individual galaxies 
indicate that there exist sub-population(s) of LBGs whose behaviors are distinct from the 
majority of LBGs in the SED-fitting. The SED-fitting generally underestimates the redshift 
slightly. 

3. The age overestimates are clearly related to the intrinsic age and specific star for- 
mation rate of each galaxy. The overestimation of mean ages is worse for galaxies with 
younger ages and higher SSFRs. Inspection of the SFHs of individual SAM galaxy confirms 
that the main origin of the bias in the age estimation is the difference of assumed SFHs in 
SAM galaxies and the simple galaxy templates used in the SED fitting. This bias in the 
agc-cstimation propagates into the stellar mass and SFR estimations, in the sense that the 
agc-ovcrcstimation leads to the mass-ovcrcstimation and SFR-undcrcstimation. The SFR- 
underestimation is further enhanced by the 'age-extinction degeneracy'. Another source of 
biases is the dominance of the current generation of star formation over the old generation(s) 
of star formation in the SED-fitting. This causes both of the stellar masses and SFRs to be 
underestimated . 

4. We perform two types of two-component SED-fitting: (1) adding a young, bursty 
component to an old component with long-lasting SFH, and (2) combining an old, bursty 
component with a younger, long-lasting component. The changes of behaviors in these 
two types of two-component fitting depend on galaxies' SFHs. Generally, compared with the 
best-fit values in the single-component fitting, the best-fit stellar masses are generally smaller 
in the two-component fitting with a young, bursty component embedded in an older, long- 
lasting component, while they become larger in the two-component fitting with an old burst 
combined with a younger, long-lasting component. For the majority of galaxies (mainly, 
with type-l/type-2 SFHs), the best-fit ages become younger and the best-fit SFRs become 
higher in both types of the two-component fitting. The behavior of galaxies with type-3 
SFH is opposite: the best-fit ages are older and the best-fit SFRs are smaller than the values 
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derived in the single-component fitting. 

5. We perform tlie SED- fitting witli different combinations of passbands - by omitting 
IRAC data or ISAAC data. If we fit tlie galaxy SEDs witli ACS and ISAAC data only 
omitting IRAC data, the derived distributions of stellar masses, SFRs and mean ages are 
significantly affected. The detailed behaviors of change in the SED-fittings with and without 
IRAC data strongly depend on galaxy's redshift and SFH. Alternatively, if we fit the SEDs 
of LBGs with ACS and IRAC data only, the derived distributions of stellar masses, SFRs, 
and means age do not show significant changes except for the enhanced bimodalities in the 
SFR/age distributions for B-/V- drop outs. This indicates that the significant changes occur 
only for a small fraction of galaxies, while the effects of ISAAC data in the SED-fitting 
are insignificant for the majority of LBGs. These experiments demonstrate the usefulness of 
the observed-frame MIR data from IRAC in constraining physical properties of the high- 2;, 
star-forming galaxies. 

6. When the allowed range of r (star formation e-folding timescale) is limited to be 
insufficiently short or long, biases in the SED-fitting increase, in general, compared with the 
case when sufficiently broad range of r is allowed (from 0.2 Gyr to 15.0 Gyr in this study). 
The mean values of the best-fit stellar mass and age are smaller than the values derived 
with the full range of r used. The age shifts are larger when we limit r to be small (< 1.0 
Gyr) than when we use very large value of r (= 15.0 Gyr). Detailed behaviors of change for 
individual galaxy depend on the SFH. 

7. The experiments with the SED templates constructed from the BC03 model (§ 15.40 
isolate the effect of limiting the allowed range of r. If r is restricted to be r < 1 Gyr, both 
the stellar mass and mean age are underestimated while the SFR is overestimated. If only 
very long values of r are used, the stellar mass is overestimated for the majority of model 
galaxies (from BC03 model) as a result of the mis-assignment of light to older stars (with 
consequently higher mass-to- light ratio). The biases in the SFR and age estimation depend 
on the age of galaxy. If the age (more exactly, value of t, i.e. time since the onset of star 
formation) is long, compatible to the age of the universe, the mean age is underestimated and 
E[B — V) is greatly overestimated leading to the severe overestimation of SFR. Otherwise, 
the mean age is overestimated and SFR is underestimated. In the case when one tries the 
single-component SED-fitting for the galaxies with two clearly distinguished generations of 
star formation (resembling repeated-burst models), the stellar mass, SFR, and mean age are 
all underestimated. 

8. Star- formation rates estimated from the UV luminosity alone may be less biased than 
those estimated from SED-fitting, provided one has a reasonable estimate of E{B — V). This 
is mainly due to the fact that the results are less subject to the degeneracy between age and 
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dust-extinction when we use only rest-frame UV photometry. However, the bias depends on 
galaxies' redshift, and more significantly on the estimation of mean dust-extinction, which is 
challenging. A relatively small change in E{B — V) would result in large bias in UV-derived 
SFRs. 

9. We show that biases arising in the SED-fitting procedure can affect studies of the high- 
redshift, star-forming galaxies. The different directions and amounts of biases depending on 
galaxy's SFH can produce the artificial bimodalities in the age or SFR distributions. This 
can affect the interpretation of the properties and nature of these galaxies. Also, the stellar 
mass underestimation for massive LBGs, combined with the SFR underestimation, can affect 
the interpretation of possible evolutionary paths for these massive LBGs. 

In conclusion, we show that single-component SED-fitting generally slightly underes- 
timates the LBGs' stellar mass distributions, while the SFR distributions are significantly 
underestimated and the age distributions are significantly overestimated. The main causes 
of these biases are: (1) the difference of assumed SFHs between in the SAM galaxies and 
simple templates used in the SED-fitting, and (2) the effects of the current generation of star 
formation masking the previous generation(s) of stellar population. The well-known 'age- 
extinction degeneracy' (or 'age-extinction-redshift degeneracy' in the case when redshifts are 
allowed to vary freely as an additional free parameter during the SED-fitting procedure) 
plays an additional role, mainly in the estimation of SFR distributions. Consequently, the 
directions and amounts of the biases in the SED-fitting strongly depend on galaxy's star 
formation history (SFH). If we change various inputs or fitting parameters in SED-fitting, 
such as the range of r - the e-folding timescale of star formation, combinations of pass- 
bands used, or assumed SFHs, the derived distributions of best-fit stellar masses, SFRs and 
ages can change dramatically. Due to the compensating causes of biases, the best-fit stellar 
mass distributions are more stable against these changes than the SFR/age distributions. 
Moreover, the behaviors of individual galaxy in various settings of the SED-fitting strongly 
depend on galaxy's SFH. 

These biases arising in the SED-fitting can have significant effects in the context of 
the galaxy formation/evolution studies as well as cosmological studies. Besides the effects 
discussed in § [3, biases in the estimation of stellar mass and SFR (which are dependent on 
the various input settings in the SED-fitting, on available passbands and more importantly 
on individual galaxy's SFH) can affect, for example, the estimation of the stellar mass 
function as well as the global density of the stellar mass and SFR, which are important 
probes in galaxy formation/evolution. The blind comparison among various works - which 
were done with different input settings in the SED-fitting, with the different sets of available 
photometric data (e.g. works done before or after Spitzer-eia), or with differently-selected 
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galaxy samples - can also misleads us. Therefore, appropriate caution should be applied to 
the estimates of physical parameters of high-redshift, star-forming galaxies through the SED- 
fitting, and also to the interpretation in the context of galaxy formation/evolution based on 
(the comparisons of) the derived results. 

We acknowledge useful discussions of the issues in this paper with our many colleagues on 
the GOODS, FIDEL, and COSMOS teams, including Casey Papovich. SL thanks Myungshin 
Im for useful suggestions applied in this work. This work is supported in part by the Spitzcr 
Space Telescope Legacy Science Program, which was provided by NASA, contract 1224666 
issued by the JPL, Caltech, under NASA contract 1407, and in part by HST program number 
GO-9822, which was provided by NASA through a grant from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, 
Incorporated, under NASA contract NAS5-26555. 
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Fig. 1. — Redshift distributions of model U-, B-, & V-dropout galaxies selected through 
Lyman break color criteria as explained in § 13.11 
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Fig. 2. — Distribution of the logarithm of stellar masses for U-dropouts (left column), B- 
dropouts (middle column), and V-dropouts (right column) when galaxy redshifts are fixed 
and all of the ACS/ ISAAC / IRAC passbands are used in the SED- fitting procedure. Figures 
in the top row are the intrinsic distributions from the semi-analytic models and figures in 
the bottom row are the distributions of best-fit values from SED-fitting. Stellar masses are 
given in Mq. 
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Fig. 3. — Distribution of the logarithm of SFRs averaged over last 100 Myr (in Mq yr^^) 
for U-dropouts (left column), B-dropouts (middle column), and V-dropouts (right column) 
when galaxy redshifts are fixed and all of the ACS /ISAAC/ IRAC passbands are used. The 
top row shows the intrinsic distributions from the semi-analytic models and the bottom row 
shows the distributions of best-fit values from SED-fitting. 
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Fig. 4. — Distribution of the stellar-mass weighted mean stellar-population ages (given in 
Gyr) for U-dropouts (left column), B-dropouts (middle column), and V-dropouts (right 
column) when galaxy redshifts are fixed and all of the ACS/ISAAC/IRAC passbands are 
used for SED fitting. Figures in the top row are the intrinsic distributions from the semi- 
analytic models and figures in the bottom row are the distributions of best-fit values from 
SED-fitting. 
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Fig. 5. — Comparison of the statistical distributions of stellar masses (left column), SFRs 
(middle column), and mean ages (right column) of U-dropouts when we use ACS^ ISAAC, 
& IRAC (the second row), when we use ACS & ISAAC only (i.e. without using IRAC 
bands, the third row), and when we use ACS & IRAC only (i.e. without ISAAC, the 
fourth row), all with the redshift held fixed. The bottom row shows the case when we vary 
redshift as a free parameter. Intrinsic distributions are shown in the top row. Stellar masses 
are given in Mq, SFRs are given in Mq yr~^, and mass- weighted mean ages are given in 
Gyr. 
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Fig. 6.— 



Same as figure 5, but for B-dropouts. 
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Fig. 7.- 



Same as figures 5 and 6, but for V-dropouts. 
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Fig. 8. — Relative errors of stellar masses, SFRs, and mean ages vs. rest-frame UV (ACS 
Veoe) and optical (IRAC 4.5/im) magnitudes, and vs. rest-framge UV (-B435 — Veoe) and 
UV-optical (^775 — m3.6 ^m) colors of U-dropout LBGs in case when ACS/ISAAC/IRAC 
passbands are used and galaxy redshifts are fixed at the input values during the SED fitting. 
One object with ASFR/ SFRsam larger than 3 is excluded in figures in the middle row for 
visual clarity. 
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Fig. 9. — Relative errors of stellar masses, SFRs, and mean ages vs. rest-framge UV (ACS 
^775) and optical (IRAC 5.8fim) magnitudes, and vs. rest-frame UV (Vgoe ~ '^775) and UV- 
optical (^775 — 777.3.6 fMm) colors of B-dropout LBGs in case when AC S / 1 SAAC / 1 RAC pass- 
bands are used and galaxy redshifts are fixed at the input values. 
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Fig. 10. — Relative errors of stellar masses, SFRs, and mean ages vs. rest-frame UV (ACS 
zg^o) and optical (IRAC 5.8;um) magnitudes, and vs. rest-frame UV (^775 — 2;85o) and UV- 
optical (2:850 — ms.g ^m) colors of V-dropout LBGs in case when ACS /ISAAC /IRAC pass- 
bands are used and galaxy redshifts are fixed. 
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Fig. 11. — Relative errors of stellar masses, SFRs, and mean ages vs. rest-frame UV (ACS 
V^oe) and optical (IRAC 4.5/xm) magnitudes, and vs. rest-framge UV (-B435 — V^eoe) and 
UV-optical (^775 — m3.6 ^m) colors of U-dropout LBGs in case when ACS /ISAAC /IRAC 
passbands are used and galaxy redshifts are allowed to vary as an additional free parameter. 
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Fig. 12. — Relative errors of stellar masses, SFRs, and mean ages vs. rest-frame UV (ACS 
^775) and optical (IRAC 5.8//m) magnitudes, and vs. rest-frame UV (Vgoe — ^775) and UV- 
optical (^775 — m3.6 fj,m) colors of B-dropout LBGs in case when AC S / 1 SAAC / 1 RAC pass- 
bands are used and galaxy redshifts are allowed to vary as an additional free parameter. 
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Fig. 13. — Relative errors of stellar masses, SFRs, and mean ages vs. rest-frame UV (ACS 
Zs5o) and optical (IRAC 5.8/im) magnitudes, and vs. rest-frame UV (2775 — ;Z85o) and UV- 
optical (-2850 — iTT-A.b ^im) colors of V-dropout LBGs in case when ACS /ISAAC /IRAC pass- 
bands are used and galaxy redshifts arc allowed to vary as an additional free parameter. 
One object with AM^/M^^sam larger than 2.5 is excluded from the figures in the top row 
for visual clarity. 
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Fig. 14. — Relative errors of stellar masses, SFRs, and mean ages vs. intrinsic stellar masses, 
SFRs, mean ages, and SSFRs for U-dropout LBGs. ACS /ISAAC /IRAC passbands are 
used and redshifts are fixed. One object with ASFR/ SFRsam larger than 3 is excluded 
from the figures in the middle row for visual clarity. 
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Fig. 15. — Relative errors of stellar masses, SFRs, and mean ages vs. intrinsic stellar masses, 
SFRs, mean ages, and SSFRs for B-dropout LBGs. ACS/ ISAAC /IRAC passbands are used 
and redshifts are fixed. 
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Fig. 16. — Relative errors of stellar masses, SFRs, and mean ages vs. intrinsic stellar masses, 
SFRs, mean ages, and SSFRs for V-dropout LBGs. ACS / ISAAC / IRAC passbands are 
used and redshifts are fixed. 
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Fig. 17. — Intrinsic age vs. specific SFR (SSFR) for B-dropout galaxies. Blue dots 
represent galaxies with 0.0 < AAge/AgesAM < 0.75. Green dots are for galaxies with 
0.75 < AAge/AgesAM < 2.0, and larger, red dots are for ones with AAge/AgesAM > 2.0. 
Purple crosses show age and SSFR for galaxies whose mean ages are underestimated, unlike 
the majority of B-dropouts. Ages are given in Gyr, and SSFRs are given in yr^^. 
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Fig. 18. — Intrinsic star formation histories of typical dropout galaxies among the galaxies 
with the largest overestimation of mean stellar population ages (blue line), along with star 
formation histories of the best-fit BC03 template (red line). These galaxies are the ones with 
type-2 SFH (see text). The red, dotted vertical line shows the point where lookback time is 
100 Myr. SFR are measured over past 100 Myr timespan in this work. Lookback time is 
given in Gyr. 
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Fig. 19. — Intrinsic star formation histories of typical dropout galaxies among the galaxies 
with the smallest overestimation of mean stellar population ages (blue line) along with star 
formation histories of the best-fit BC03 template (red line). These galaxies are the ones with 
type-1 SFH (see text). The red, dotted vertical line shows the point where lookback time is 
100 Myr. SFR are measured over past 100 Myr timespan in this work. Lookback time is 
given in Gyr 
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Fig. 20. — Correlations among relative errors of stellar masses, SFRs, and mean ages for 
U-,B-,V-dropout galaxies. AC S / 1 SAAC / 1 RAC passbands are used and redshifts are fixed 
in the SED-fitting. 
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Fig. 21. — Relative errors of stellar masses, SFRs, and mean ages vs. intrinsic stellar masses, 
SFRs, mean ages, and SSFRs for U-dropout LBGs. ACS/ISAAC/IRAC passbands are 
used and redshifts are allowed to vary as a free parameter. 




Fig. 22. — Relative errors of stellar masses, SFRs, and mean ages vs. intrinsic stellar masses, 
SFRs, mean ages, and SSFRs for B-dropout LBGs. ACS/ ISAAC /IRAC passbands are used 
and redshifts are allowed to vary. 
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Fig. 23. — Relative errors of stellar masses, SFRs, and mean ages vs. intrinsic stellar masses, 
SFRs, mean ages, and SSFRs for V-dropout LBGs. ACS/ISAAC/IRAC passbands are 
used and rcdshifts arc allowed to vary. One object with AM^/M^^sam larger than 2.5 is 
excluded from the figures in the top row for visual clarity. 
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Fig. 24. — Correlations among relative errors of stellar masses, SFRs, and mean ages for U- 
,B-,V-dropout galaxies. ACS/ISAAC/IRAC passbands are used and redshifts are allowed 
to vary. 
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Fig. 25. — Star formation histories of typical B-dropout galaxies with underestimated mean 
stellar population ages (type-3 SFH, see text). The dotted vertical hne shows the point 
where lookback time is 100 Myr. SFR are measured over past 100 Myr timespan in this 
work. Lookback time is given in Gyr. 
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Fig. 26. — Stellar mass distributions of U-dropouts (left column), B-dropouts (middle col- 
umn), and V-dropouts (right column). Figures in the top row show the distributions of 
intrinsic stellar masses. The middle row is for stellar masses from the two-component fitting 
with a secondary young burst added to a main component (as investigated in § IS.l.ip . and 
the bottom row is for stellar masses from the two-component fitting with a maximally old, 
burst component plus a younger, long-lasting component (§ I5.1.2p . Stellar masses are given 
in Mm. 
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Fig. 27. — SFR distributions of U- drop outs (left column), B-dropouts (middle column), and 
V-dropouts (right column). Figures in the top row show intrinsic SFR distributions. The 
middle row is for SFRs from the two-component fitting with a secondary young burst added to 
a main component (§ l5.1.ip . and the bottom row is for SFRs from the two-component fitting 
with a maximally old burst component plus a younger, long-lasting component (§ I5.1.2p . 
SFRs are given in Mq yr~^. 
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Fig. 28. — Mass- weighted stellar-population mean age distributions of U-dropouts (left col- 
umn), B-dropouts (middle column), and V-dropouts (right column). Figures in the top row 
show intrinsic age distributions. The middle row is for ages from the two-component fitting 
with secondary a young burst added to a main component (§ 15.1.11) . and the bottom row 
is for ages from the two-component fitting with a maximally old burst component plus a 
younger, long-lasting component (§ [5lL2])- Age is given in Gyr. 
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Fig. 29. — Correlations between relative errors in two-component fitting with a secondary 
young burst (§ IS.l.ip and relative errors in single- component fitting. Figures in the top 
row show correlations between relative age errors in two-component fitting and relative age 
errors in the single-component fitting. The middle row is for correlations between relative 
stellar mass errors in the two-component fitting and relative stellar mass errors in the single- 
component fitting. The bottom row is for correlations between relative SFR errors in the two- 
component fitting and relative age errors in the single-component fitting. A"'value represents 
^ (value from two-component fitting) - (intrinsic value)\ and A^value represents ^ (value from 
single- component fitting) - (intrinsic value)\ 
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Fig. 30. — Ratio of the best-fit stellar masses ((a)-(c)), SFRs ((d)-(f)), and mean ages ((g)- 
(i)) with and without IRAC photometry for U-dropouts ((a), (d), and (g)), B-dropouts ((b), 
(e), and (h)), and V-dropouts ((c), (f), and (i)) dependent on the best- fit values with IRAC 
photometry. 'AS' stands for values without IRAC photometry, and 'ASR' for values with 
IRAC photometry. Objects with AgtAsl-^g^ASR > 4.5 are excluded from the figures in the 
bottom row for visual clarity. 
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Fig. 31. — {AgCri — AgCaii) / AgCaii as a function of intrinsic age (top row), intrinsic specific 
SFR (middle row), and {Agcaii — AgesAAi) / AgesAM (bottom row) for U-dropouts (left col- 
umn), B-dropouts (middle column), and V-dropouts (right column). Here, Agcri and Agcaii 
are the best- fit mean ages derived when we limit r as < 1.0 Gyr and when we allow r to 
vary from 0.2 Gyr to 15.0 Gyr, respectively. AgesAM is an intrinsic age, and SSFRsam is 
an intrinsic SSFR. 
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Fig. 32. — Best-fit tau dependent discrepancies in best-fit t (t^-is — taii). taii is the best-fit t 
obtained when we allow the full range of r, and t^-is is the best-fit t derived if we use single 
value of T (= 15.0 Gyr) in the SED-fitting. ^,-15 and taii are in Gyr. 
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Fig. 33. — Ratio of SFR derived from rest-frame UV luminosity to intrinsic SFR as a function 
of redshift for B-dropouts. All B-dropouts are assumed to be at 2; = 4.0, and to be extincted 
by dust with amount of E{B — V) — 0.15. 
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Fig. 34. — Distributions of intrinsic SFRs {left), and SFRs derived from rest-frame UV 
luminosity assuming all galaxies are aX z = 4.0 and have mean dust-extinction of E{B — V) — 
0.15 (right) for B-dropout galaxies. 
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Table 1. Fitting Parameters 



IMF 


T (Gyr) t (Gyr) metallicity {Zq) 


internal extinction 


IGM extinction 


Chabrier 


0.2-15.0 0.01-2.3^ 0.2, 0.4, 1.0 


Calzetti'' 


Madau 



H is limited to be smaller than the age of the universe at each galaxy's redshift. 
^0.000 < E{B -V) < 0.950 with AE{B -V) = 0.025 
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Table 2. Mean Values of Stellar Population Parameter Distributions 





Rprl<jhif(- 

JTLtlLlbllliL 


OttliiCll iV±Clk5& I iK/0 J 






U-dropouts 










SAMs^ 


3.36 


9.376x109 


11.005 


0.465 


SEDz-fix^ 


3.36 


7.594x10^ 


3.887 


0.915 


S ED z- free' 


3.27 


7.038x10'^ 


4.682 


0.870 


B-dropouts 










SAMs^ 


4.02 


8.888x10^ 


15.650 


0.353 


SEDz-fix^ 


4.02 


6.680x10^ 


6.638 


0.667 


SEDz-free^ 


3.88 


4.327x10'' 


19.891 


0.420 


V-dropouts 










SAMs*^ 


4.97 


7.946x10^ 


18.920 


0.249 


SEDz-fix^ 


4.97 


5.983x10^ 


7.287 


0.508 


SEDz-free^ 


4.88 


4.496x10^^ 


16.905 


0.344 



^Intrinsic values from SAM catalogs 

'^Values derived in the SED-fitting with redshifts fixed as input values 
■^Values derived in the SED-fitting with redshifts allowed to vary freely 
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Table 3. Means and Standard Deviations of Relative Offsets^ for various Stellar 

Population Parameters 



Redshift Stellar Mass SFR Mean Age 



U-dropouts 

SEDz-free^ 

B-dropouts 
SED, 
SED, 



-fix 



-free 

V-dropouts 



SED 



z—fix 



SED 



-0.021±0.015 



-0.030±0.031 



-0.039±0.232 
-0.108±0.230 



-0.120±0.228 
-0.371±0.355 



-0.585±0.185 
-0.530±0.258 



-0.558±0.288 
0.212±1.181 



1.075±0.510 
0.955±0.502 



1.009±0.586 
0.273±1.012 



z—free 



■■■ -0.169±0.179 -0.597±0.133 1.162±0.589 
-0.014±0.012 -0.339±0.274 -0.094±0.964 0.487±0.875 



^Relative offset is defined as (V aluesED — V aluesAKi) / iV dluesAivi) for 
stellar mass, SFR, and age. For redshift, relative offset is defined as {zsed— 
zsam)/{^+zsam)- 

^Relative offsets between values derived in SED-fitting with redshifts 
fixed and intrinsic values 

■^Relative offsets between values derived in SED-fitting with redshifts 
allowed to vary and intrinsic values 
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Table 4. Relative Changes'^ of the Mean Stellar Masses and Ages with the Limited r 

Range 



rs Used for Fitting Stellar Mass (%) Mean Age (%) 



U-dropouts 

r < 1.0 Gyr -4.9 -15 

r = 15.0 Gyr -8.0 -12 

B-dropouts 

r < 1.0 Gyr -12 -26 

r = 15.0 Gyr -3.4 -7.4 

V-dropouts 

T < 1.0 Gyr -5.4 -16 

T = 15.0 Gyr -3.2 -6.6 



Relative change is defined as ({valuer) ~ 
(valueaii))/ {value all). {valueau) is the mean value 
of the stellar mass or age derived with the full range of 
T. {valuer) is the mean value of the stellar mass or age 
derived with the limited r range. 
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Table 5. Toy Models 



Model 1 Model 2 Model 3 



input 

T (Gyr) 15.0 0.2 0.2 (young), 0.2 (old) 

t (Gyr) 0.01, 0.05, 0.1, 0.5, 1.0 0.1, 0.2, 0.5, 1.0, 1.3 0.1 (young), 1.0 (old) 

E{B - V) 0, 0.05, 0.1, 0.2 0, 0.05, 0.1, 0.2 0.0, 0.2 

Mass Fraction^ ■ ■ ■ ■ ■ ■ 0.1, 0.26, 0.52, 1.0, 2.6 

Number of SEDs 120 120 60 

used for fitting 

r (Gyr) < 1.0 > 8.0 0.2 < r < 15.0 

t (Gyr) 0.01 - 2.3*^ 0.01 - 2.3*^ 0.01 - 2.3*^ 

'^Stellar mass fraction between young component and old component defined as 

M young/ Mold 

H is fimited to be smaller than the age of the universe at each galaxy's redshift. 



